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- - ■ 1. Introduction 

The structure of the nucleon in the non-perturbative regime of QCD can be characterized 
by a collection of scales, related to the probe which is used to study this basic building 
block of atomic nuclei. The typical nucleon size as deduced from electron scattering 
experiments is about 0.85 fm, whereas model-dependent analyses of proton-antiproton 
annihilation lead to a baryon charge radius of about 0.5 fm. An intermediate scale 
is set by the weak charged currents, which allows to measure the axial charge radius, 
rA — 0.65 fm, and related moments. For a more detailed discussion of these various 
scales and how they can be understood in a simple model, see reference [Q. Here, we will 
be concerned with nucleon matrix elements of the axial current, giving a status report 
on our current knowledge about the pertinent coupling constants and form factors, 
with particular emphasis on the theoretical methods which are used to extract this 
information from experiment. Our basic object is the QCD axial current, expressed in 
terms of the light quark fields 

A; = q^,^,T'^q, (1) 
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where for the two and three flavor case, the pertinent quark fields, generators and flavor 
indices are given by 
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with T°' the conventional Pauli isospin matrices and the A" are Gell-Mann's SU(3) 
matrices. Note that we do not consider the singlet components corresponding to a = 
in what follows. To be specific, consider the matrix-element of the SU(2) isovector axial 
quark current between nucleon states. 



{N{p') |A^ \N{p)) = u{p') 



l,GA{t) + 



jp' -P)fi 
2m 



Gp{t) 



l^^u{p) , 



(3) 



with t = {p' — pY the invariant momentum transfer squared and m = {rrip + mn)/2 
the nucleon mass. Note that in this review, we mostly consider the isospin symmetric 
case and only differentiate between the proton {rrip) and the neutron (m„) masses in 
kinematical factors. The form of equation (^) follows from Lorentz invariance, isospin 
conservation and the discrete symmetries C, P and T and the absence of second class 
currents 0, which is consistent with experimental information, see e.g. reference 0]. 
Gyi(t) is called the nucleon axial form factor and Gp{t) the induced pseudoscalar form 
factor. Consider first the axial form factor, which probes the spin-isospin distribution 
of the nucleon (since in a non-relativistic language, this is nothing but the matrix- 
element of the Gamov-Teller operator cr r). The axial form factor admits the following 
expansion for small momentum transfer 



GAit) = QA 



i + g(ri)t + o(t^; 



(4) 



with qa the axial-vector coupling constant, qa = 1.2673 ± 0.0035 0, and (r^)^''^ the 
nucleon axial radius. The experimental determinations and theoretical understanding 
of this fundamental low-energy observable will constitute a main part of this review. 
The axial-vector coupling constant is measured in (polarized) neutron beta-decay, for 
recent references see e.g. [^, Q. Here, we will not discuss the determinations of qa and 
of the related CKM matrix-element Vud, for a review see [0]. For low and moderate 
momentum transfer, say |t| < 1 GeV^, the axial form factor can be represented by a 
dipole fit 

9A 



GAit) 



t/Mj^' ' 



(5) 



in terms of one adjustable parameter. Ma, the so-called axial mass (or sometimes dipole 
mass). Therefore, one can express the axial radius in terms of the dipole mass 
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The induced pseudoscalar form factor is believed to be understood in terms of pion pole 
dominance and we will discuss the theoretical status concerning the corrections to this 
leading order term. Most measurement of this form factor stem from ordinary muon 
capture on the proton, /i~ + p — i> z/^ + n. More precisely, in that reaction one measures 
the induced pseudoscalar coupling constant gp, 

9p = ^Gpit = -0.88M'^), (7) 

with M^ the muon mass. This coupling constant is nothing but the value of the induced 
pseudoscalar form factor Gp{t) at the four-momentum transfer for muon capture by 
the proton at rest. 



t = -M^ 



1- 



(M^ + THp) 



2-m2 



,w.. N -0.88 M^. (8) 

In what follows, we will review these determinations and also discuss radiative muon 
capture as well as pion electroproduction which allow to investigate the momentum 
dependence of Gp{t). 

The material is organized in the following way. In section ^, we collect the available data 
on the axial form factor Ga(^), the induced pseudoscalar coupling constant gp and the 
corresponding form factor Gp{t). Section]^ is devoted to the theoretical methods which 
allow to determine the axial form factor and the axial radius of the nucleon as well as 
the pion charge radius from (anti)neutrino scattering and pion electroproduction data. 
The QCD analysis of the induced pseudoscalar form factor and the methods to extract 
it from ordinary and radiative muon capture are discussed in section ^ Section ^ is 
devoted to a discussion of certain axial matrix elements in the baryon octet pertinent to 
the quark mass analysis. A short summary and outlook are given in section ^. Before 
proceeding, we point out that in this review no model calculations of the matrix elements 
are reviewed, as illustrative as they might be. Our emphasis is on applying a model- 
independent approach tightly constrained by the symmetries of the Standard Model to 
determine the axial structure of the nucleon in the non-perturbative regime. 

2. Data overview 

In this section, we will give a brief review of published data on the axial and induced 
pseudoscalar form factors. The methods which lead to these results will be presented 
and critically discussed in the following sections. At this point, we only wish to point 
out that systematical errors for certain processes have certainly been underestimated, 
but here we simply collect the data as they are available. 

2.1. Data on the axial form factor 

There are two methods to determine the axial form factor of the nucleon, namely 
(anti)neutrino scattering off protons or nuclei and charged pion electroproduction. 
We consider first the experimental data coming from measurements of (quasi)elastic 
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anti)neutrino scattering off protons p, §, |T^, off deuterons |TT[]-|T^ and otlier nuclei (Al 



Fe) [0, [T8| or composite targets like freon [^-[^ and propane |2^, |23[. In the left panel 
of figure [l| we show the available values for the axial mass Ma obtained from neutrino 
scattering experiments. As pointed out in [^], references Jl^ |T^, ^, ^ reported 
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Figure 1. Axial mass Ma extractions. Left panel: From (quasi) elastic neutrino 
and antineutrino scattering experiments. The weighted average is AIa = (1-026 ± 
0.021) GeV. Right panel: From charged pion electroproduction experiments. The 
weighted average is Ma = (1.069 ± 0.016) GeV. Note that value for the MAMI 
experiment contains both the statistical and systematical uncertainty; for other values 
the systematical errors were not explicitly given. The labels SP, DR, FPV and BNR 
refer to different methods evaluating the corrections beyond the soft pion limit as 
explained in the text. 



severe uncertainties in either knowledge of the incident neutrino flux or reliability of the 
theoretical input needed to subtract the background from genuine elastic events (both 
of which gradually improved in subsequent experiments). The values derived in these 
papers fall well outside the most probable range of values known today and exhibit 
very large statistical and systematical errors. Following the data selection criteria of 
the Particle Data Group [^, they were excluded from this compilation. In all cases, 
the axial form factor data were parameterized in terms of a dipole, the resulting world 



average is 



Ma = (1.026 ± 0.021) GeV (neutrino scattering) . 



(9) 



The other determinations of the axial form factor are based on the analysis of charged 
pion electroproduction off protons, see references ||2^[ill-[lll, slightly above the pion 
production threshold (note that the MAMI measurement is presently extended p5| to 
lower momentum transfer and to check the cross section at the highest Q^ point reported 
24|). Such type of analysis is more involved. It starts from the low-energy theorem of 



m 



Nambu, Lurie and Shrauner 



37| for the electric dipole amplitude Eq^' at threshold. 
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valid for soft pions, i.e. pions with vanishing four-momentum. Model-dependent 

corrections (so-called hard pion corrections) were developed to connect the low-energy 

theorem to the data, that is to the real world with a finite pion mass, see references [pH] - 

|3], labeled SP, FPV, DR and BNR, respectively. For a given model, the values of 



the axial mass were determined from the slopes of the angle-integrated differential 
electroproduction cross sections at threshold. The results of various measurements and 
theoretical approaches are shown in the right panel of figure |l]. Note again that references 



2% |38| were omitted from the fit for lack of reasonable compatibility with the other 
results. In figure § we have collected the various electroproduction data, in comparison 
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Figure 2. Experimental data for the normalized axial form factor extracted 
from pion electroproduction experiments in the threshold region. Note that for 
the experiments where various theoretical models were used in the analysis to 
extract Ga, all results are shown. For orientation, the dashed line shows a dipole 
fit with an axial mass Ma = 1.1 GeV. 



to a dipole fit with Ma = 1.1 GeV. Again, at various values of the momentum transfer 
one sees two or three data points, these show the model-dependence due to the applied 
hard pions corrections. The resulting world average of the dipole masses collected in 
figure is 

Ma = (1.069 ± 0.016) GeV (electroproduction) . (10) 

Although some of these results have large uncertainties, the weighted average for the 
axial mass determinations from neutrino scattering and pion electroproduction are quite 
precise. In particular, we notice an axial mass discrepancy, i.e. the so determined axial 
masses differ significantly, 

AMa = (0.043 ± 0.026) GeV , (11) 
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which translates into an axial radius difference of about 5%, i.e. the axial radius of 
the nucleon determined from neutrino scattering data appears bigger than the one from 
electroproduction, {r\)l^^,^^,^^ = (0.666±0.014) fm versus {r\)l(^^^^ = (0.639±0.010) fm, 
the axial radius discrepancy to be discussed in section H. 



2.2. Data on the induced pseudoscalar form factor 

Ordinary muon capture (OMC), 

/.-(/) +p(r)^z/^(0+n(r'), (12) 

where we have indicated the four-momenta of the various particles, allows to measure the 
induced pseudoscalar coupling constant, gp. Most experiments so far have used liquid 
hydrogen targets, which leads to atomic physics complications because the /xp atom 
rapidly forms an p/ip orthomolecule. Consequently, one has to know the transition 
rate from the ortho to the para ground state. Furthermore, there can be mixing 
between S = 1/2 and S = 3/2 components of the ortho p//p state fi^. To facilitate 
the discussion, the various atomic and molecular states are shown in figure |^ together 
with the pertinent transition rates (as determined experimentally). The world data on 
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Figure 3. Level scheme for OMC in liquid hydrogen. In the atomic (p/i) state, 
one has a triplet/singlet state where the sum of the proton and the muon spins is 1/0 
(upper two levels). In the ortho (pup) state, the two proton spins sum up to one, 
splitting into two levels with S = 3/2 and S = 1/2 as indicated. Weinberg's mixing 
suggestion Q applies to these two levels. 



electronics experiments for muon capture in hydrogen [0-[0 are collected in figure ^, 
these values have been rescaled to the present day values of the axial-vector coupling 
constant, the Fermi constant and the Cabbibo angle. The weighted world average is 

gp = 8.79 ± 1.92 , (13) 

which is in good agreement with theoretical expectations as discussed below. At this 
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Figure 4. Determination of the pseudoscalar coupling constant gp from (electronics) 
experiments. The numbers have been rescaled to the present value of gA = 1.267. For 
the measurements on liquid hydrogen, the rate Aop = (4.1 ± 1.4) • 10''' s~^ is assumed. 
The weighted average is gp — 8.79 ± 1.92. 



point it is worth to emphasize that the experimental value for the ortho-para transition 
rate Aqp = (4.1 ± 1.4) ■ 10^ s^^ |^ is used to arrive at the world average. This 
value disagrees, however, with the theoretical expectation of reference ||5^, Aqp = 
(7.1 ± 1.2) • 10^ s^^. We also point out that the theoretical formula for the liquid capture 



rate a la |5^ involves a time integration from zero to infinity while in the measurement 
only a finite time interval is considered, thus a direct comparison is dangerous. A similar 



observation has recently been made in reference |5J]. Stated differently, the theory has 
to adopt to the experimental circumstances. We would like to stress here the importance 
of the new proposed experiment at PSI |55| which will be done with a hydrogen gas 
target and will thus be independent of possible molecular complications because one 
will directly measure the atomic singlet rate. This experiment attempts to measure 
gp to one percent accuracy, which would be a significant improvement compared to 
the present situation. These topics will be taken up again in section |. The induced 
pseudoscalar coupling constant has also been determined by muon capture on various 



nuclei, from light to heavy ones. A comprehensive review has been given in [Q, here 
we only report results for some light nuclei in table |1]. Note that these determinations 
are mostly in agreement with the theoretical expectations, but the spread in the results 
is larger than given in the table ||56[. In heavier nuclei, there has been much discussion 
about a possible quenching of gp (similar to the quenching of gA), but the situation 
can not be considered settled. We do not pursue this topic here, but refer to the recent 
reference [^. While the momentum transfer in OMC is fixed, radiative muon capture 
(RMC), 

fj.-{l) + p{r) -^ v^{l') + n{r') + ^{k) , (14) 
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Table 1. Pseudoscalar coupling constant determined from OMC in light nuclei. 
Nucleus gp Reference 



^He (capture to triton) 8.6 ± 1.5 

^^C (capture to ground state) 8.3 ± 2.5 
1^0 (capture to isN(0")) 10.0 ± 1.2 




has a variable momentum transfer t and one can get up to t = M^ at the maximum 
photon energy of about k^ ~ 100 MeV, which is quite close to the pion pole. This 
amounts approximately to a four times larger sensitivity to gp in RMC than OMC. 
However, this increased sensitivity is upset by the very small partial branching ratio in 
hydrogen (~ 10~^ for photons with fco > 60 MeV) and one thus has to deal with large 
backgrounds. Precisely for this reason only very recently a first measurement of RMC 
on the proton has been published [^ ^. The resulting number for gp, which was 
obtained using a relativistic tree model including the A-isobar ||6^ to fit the measured 



photon spectrum, came out significantly larger than expected from OMC, 

^RMC ^ ^2.35 ± 0.88 ± 0.38 ~ 1.4 g^^^ , (15) 

and thus also about 40% above all theoretical expectations (see section [4.1| ). It should 
be noted that in this model the momentum dependence in Gp{t) is solely given in terms 
of the pion pole and the induced pseudoscalar coupling is obtained as a multiplicative 
factor from direct comparison to the photon spectrum and the partial RMC branching 
ratio (for photon energies larger than 60 MeV). We will critically examine this procedure 
in section [4.3.2| . It was also argued in |^ |6^ that the atomic and molecular physics 
related to the binding of the muon in singlet and triplet atomic fip and ortho and 
para pfip molecular states is sufficiently well under control (which is under debate, see 
section |i7^ ). This result spurred a lot of theoretical activity, as discussed later, but 
here we only wish to point that it is a viable possibility that the discrepancy does not 
come from the strong interactions but rather is related to the distribution of the various 
spin states of the muonic atoms. It is therefore mandatory to sharpen the theoretical 
predictions for the strong as well as the non-strong physics entering the experimental 
analysis. In principle, Gp{t) can also be measured in pion electroproduction, because it 
enters the longitudinal cross section (in parallel kinematics) together with many other 
effects (see the discussion in section p.4.4|) . However, the leading dependence due to the 
pion pole (and its chiral corrections) is unique and can be tested. So far there has only 
been one experiment |^ that took up the challenge, with the results shown in figure |. 
The observed momentum dependence agrees with theoretical expectations, but a more 
refined measurement is certainly called for. In fact, at the Mainz Microtron MAMI-B 
a dedicated experiment has been proposed to measure the induced pseudoscalar form 
factor by means of charged pion electroproduction at low momentum transfer [Q. At 
this point, it is important to stress that in the low-energy region, Gp{t) is not truly 



independent from Gyi(t), as discussed in detail in section ETI. We conclude that the 
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induced pseudoscalar form factor is the least well known of all six electroweak nucleon 
form factors. 
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Figure 5. The "world data" for the induced pseudoscalar form factor Gp{Q^). 
The pion electroproduction data (filled circles) are from reference |p5[. Also shown 
is the world average for ordinary muon capture at Q^ — 0.88A/^ (diamond). For 
orientation, we also show the theoretical predictions discussed later. Dashed curve: 
Pion-pole (current algebra) prediction. Solid curve: Next-to-leading order chiral 
perturbation theory prediction. 



3. Nucleon axial radius and form factor 

3.1. Determination of the axial form factor I: Neutrino scattering formalism 

Neutral and charged current (anti)neutrino scattering off protons or light nuclei can 
be used to extract the axial form factor. Since these topics are well documented in 
the literature [16 



d8|, ^, we only give some basic formulae (sometimes in a highly 
symbolic notation). The starting point is the effective Lagrangian of the electroweak 
Standard Model for elastic neutrino-hadron scattering, such as u^p -^ u^p, or quasi- 
elastic scattering, such as z/^n — > fi^p reactions (and similarly for anti-neutrinos), 

Grp 



^uh ^ _r^ 



^£7'^(1+75)^J^ 



hadr 



+ h.c. 



(16) 



with I the neutrino/lepton in the final state. We are interested here in the hadronic 



current J]]^^ which can be written as 



J. 



hadr 



aV'^+pAl + ^V'^+5Al + ... 



(17) 



where the up- and down-quark contributions are combined to form the isoscalar and 
isovector vector (V^, V^) and axial-vector currents (y4°,y4^). The ellipses represent 
heavy quark (s, c, h) terms. In the Standard Model, one has at tree level 
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1 - 2 sin^ Ow , /? = 1 
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with 6]Y the (weak) mixing angle. It describes the strength of the mixing of the 
electromagnetic current with the neutral weak current and thus only appears in the 
vector currents. Assuming now time reversal invariance, isospin invariance, the absence 
of second class currents, and neglecting the induced pseudoscalar form factor (in 
some analyses, this form factor is retained in its simplified pole form, Gp{Q'^) ~ 
2mG A^Q"^) / iQ'^ + M^)), the hadronic current matrix element between nucleon states 
is given by 



\N{p)) = u{p') 



i,F,m + 



2m 



F2{Q')+1,15Ga{Q' 



u{p) , (19) 



-q^q'^ the squared momentum transfer. Because of vector 



with Q2 = 

current conservation, the Dirac and Pauli vector form factors appearing in equation (p!9|) 
are the same as the ones measured in elastic electron-hadron scattering and can 
therefore be considered as known. What is usually assumed in the analysis of the 
neutrino scattering experiments is the dipole parameterization for the Sachs form factors, 
Ge = Fi + (QV4m2)F2 and Cm = Fi + Fg. It has the form: G^ = D, G^ = 0, 
^M — f^p/nD, where D = {1 + Q'^/My)~'^, Kp/n is the proton/neutron magnetic 
moment and My = 0.84 GeV. The axial form factor is then extracted by fitting the 
Q^-dependence of the (anti)neutrino-nucleon cross section. 



dQ' 



SnE,, 



AiQ')TBiQ' 



[s -u 



+ GiQ' 



u 



(20) 



771^ m^ 

in which Gyi(Q^) is contained in the bilinear forms A{Q'^), B{Q'^) and C{Q'^) of the 
relevant form factors and is assumed to be the only unknown quantity 



A 
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4Q2 
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'A 



;i - r)Fl + r(l 



r)F| 



M, 



AtF.F, - -f ((Fi + F^y + G\) 



r)G' 
m^ [■ . . . ^ . . - ^ 

AtGa{Fi + F2) , 

C = \[Gl + F^ + rF^] , 

where the plus sign is for neutrinos and the minus sign for antineutrinos, and 



(21) 



u 



AmEy -Q , T 
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with Q^ = Q'^ — M't and M^ the muon mass. Of course, for elastic scattering u p 



(22) 

(-) 
► u p 

the terms ~ M^ should be dropped. Due to the smallness of the lepton mass, they are 

usually also dropped for the quasi-elastic reactions. Furthermore, E^, is the neutrino 

energy which is usually determined from the neutrino flux spectrum. Using the the 

dipole parameterization equation (^) allows then to determine Ma from a flt to the 

measured (anti)neutrino cross section. Here, one assumes that heavy quark corrections 

to the axial current can be neglected, although they are not necessarily small ||7D|, [71 



In particular, the so-called strange quark content of the nucleon has attracted much 
attention over the last decade. Since the Q^-dependence for such type of contribution 
is not known, one can simply modify the dipole ansatz to Ga{Q^) = Gdipoie(Q^)(l + ^), 
with rj parameterizing the strange (heavy) quark contribution. We do not wish to further 



elaborate on this topic but rather refer to the references ||10|, [72 
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3.2. Determination of the axial form factor II: Electroproduction formalism 
Consider now pion electroproduction off nucleons, 

7^(^:2) + iVi(pi) ^ n'^iq) + N2{p2) , (23) 

where 7* denotes the virtual photon with virtuality k"^ < 0, Ni{pi) {i = 1,2) the 
initial/final nucleon and vr"(g) the pion with Cartesian isospin a = (0, +, — ) and four- 
momentum q^. We will also use the positive definite quantity Q^ = —k"^. The pertinent 
Mandelstam variables are s = (pi + fc)^, t = {pi — P2Y and u = {pi — qY, subject 
to the constraint s + t + u = 2m^ + M^ + k"^, with M^r the pion mass. In the Born 
(one-photon-exchange) approximation, the corresponding coincidence cross section can 
be factorized as |7^ 

± = Y ^ (2A) 

dE'^d^'.dni ^dn%' ^ ' 

where T^ is the virtual photon flux, E'^, Q'^ the energy and the solid angle of the scattered 
electron, and da^/dQ* is the virtual photon cross section in the centre-of-mass frame 
of the final ttN system, as denoted by the star. It can be further decomposed into 
transverse, longitudinal and two interference parts, 

da^ da-T ^ doi^ /n *n , ^ '^^lt , , da^T „ , ,^^. 

d^=rff]^ + ^^rf^ + V2^i^(l + ^)7^^-^'^'^ + ^^^-^2</>. (25) 

with the transverse (e) and longitudinal (cl = —k'^e/k^) polarizations of the virtual 
photon fixed by the electron kinematics. In parallel kinematics (^* = ^^r = 0°, with ^^ 
the polar angle in the scattering plane as seen in the laboratory system) the interference 
parts vanish. Therefore, at constant four-momentum transfer, the transverse and the 
longitudinal cross sections can be separated using the Rosenbluth method by varying e, 

da^ _ daT _ ^ da^ ,^„. 

dQ* ~ dQ* ^ kf dQ*' ^ ' 

(Note that often the photon energy is denoted by u, however, here this symbol is entirely 
reserved for the pion energy). At low energies, the connection to theory is most easily 
made by means of the multipole expansion. For doing that, one considers the transition 
current related to equation (|23|), which can be decomposed in terms of six invariant 
amplitudes (we follow the conventions and notations of reference [0], see also ||75[| ) 



e • T{p2, S2; q, a\pi, Si; /c) = « M2 75 J^ e ■ -^i Ais, u) Ui , (27) 

with Si the spin index of nucleon i. The explicit forms of the operators M.i can be found 
in 



|74| and the Ai are invariant functions that depend on two kinematical variables. 
Here, we have chosen the Mandelstam variables s and u. The amplitudes Ai{s,u) have 
the isospin decomposition 

A,{s,u) = A^;-\s,u)5a, + A^r\s,u)\[T,,T,]+ Af\s,u)T, . (28) 
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In terms of the isospin components, the physical channels are given by 

r,(7'^p^A)= tw + tw, 

T,(yn^7rM= TW-T(°), (29) 

where we have listed the neutral pion channels for completeness. Most of the following 
discussions will be focused on the first reaction. For the discussion of the low-energy 
theorems that link the axial form factor to the charged pion production amplitudes, 
we spell out the corresponding multipole decomposition of the transition current 
matrix element at threshold. In the 7*N centre-of-mass system at threshold, i.e. for 
q^ = (M^, 0, 0, 0), one can express the current matrix element in terms of the two S-wave 
multipole amplitudes, called Eq^ and Lo+, 

e T = 4m{l + /x) x\ {^o+(/i, i^) + [Lo+it^, i^) - Eo+ifi, u)] k ■ ek ■ a^ Xi , (30) 

with Xij two-component Pauli-spinors for the nucleon. We have introduced the 
dimensionless quantities 

/i= , 1^ = ^- (31) 

The multipole i?o+ characterizes the transverse and Lo+ the longitudinal coupling of the 
virtual photon to the nucleon spin. For an explicit expression of these multipoles in terms 
of the Ai evaluated at threshold, Sthr = wi^(H-/i)^ andwthr = m^(l— yU— /i^+/iz/)/(H-/i), 
we refer to |7^ . As it will become clear in the following paragraphs, the two parameters 



defined in equation (|3T| ) will play the role of small expansion parameters in the QCD 
analysis of pion electroproduction. 

3.3. Current algebra derivation 

While the methods of current algebra (CA) should be considered outdated by now, we 
consider it nevertheless instructive to briefly review the CA argument which gives the 
link between charged pion electroproduction and the nucleon axial form factor. Using 
standard LSZ formalism, the transition matrix element can be written as 

T; = - j d\e"^-^{n + Ml){N2\rTx\x)V,mNi) , (32) 

with Vi, the electromagnetic current operator to which the nucleon couples, q^ the pion 
four-momentum and T is the conventional time-ordering operator. Inserting now the 
PCAC relation 9'^A°(x) = M^F^7r°(x) and performing a partial integration (note that 
while the meaning of the PCAC relation has been understood since long |]76|, ^, it does 
not offer a systematic way of calculating higher order corrections), one gets 

M'^ - q' 
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C;; = fd^xe^'^-^Sixo) {N2\[A'',{x),VM]\Ni) ■ (33) 



Using now the CA relation [^0(^)5 ^^^(0)] = ie^^^5^{x)A'l{0), one obtains in the soft pion 
hmit q^ ^ 

T: = t-- {N2{p2)\Al{0)\N,{p,)) + pole terms . (34) 

The pertinent nucleon matrix element appearing here was already given in equation (^, 
i.e. it contains the axial form factor GA{k'^) in terms of the photon virtuality A;^ < 0, 
since in the soft-pion limit q^t = the momentum transfer is {p2—pi)^i = k^. Because of 
the e-tensor, the axial form factor can only appear in charged pion electroproduction. 
In the chiral limit, Nambu, Lurie and Shrauner |36, 37] have put this in the form of a 



low-energy theorem for the electric dipole amplitude -Eo+ , 



^0+ 



a:>(M. . 0, .^) ^ ^1 - ^ 5^ {G.ie) + jJ^g'„(*^)} . (35) 



with G\.j^ the nucleon isovector magnetic form factor. The second term in this equation 
stems from the nucleon pole contribution in the chiral limit. Expanding this to 0(A;^), 
one finds 

i^«(M. = 0,fc^) = ^{l + ^(r^) + |!-(K. + i)+O(A:3)}, (36) 

with hiy the nucleon isovector anomalous magnetic moment. Of course, this equation is 
only correct for massless pions with zero three-momentum. Before the advent of chiral 
perturbation theory, there existed many prescriptions to get a handle on the corrections 



due to the pion mass, see references [^-[^ For a textbook treatment of this issue, see 



e.g. reference [73|. This was already addressed in section ^ in connection with theoretical 
uncertainties of the axial form factor extracted from electroproduction based on various 
prescriptions to deal with the corrections beyond the soft-pion limit. These uncertainties 
can and have to be overcome using the modern methods discussed next. To end this 
section, we remark that in principle equation (0) should also include the induced 
pseudoscalar form factor. However, in soft-pion limit the Gp(t) term is proportional to 
the photon momentum A;^, and when contracted with the photon polarization vector, 
one obtains the structure e ■ kGp(t). Therefore, the induced pseudoscalar form factor 
only contributes to the longitudinal cross section, more precisely to the longitudinal 



S-wave multipole Lo+ (as discussed in section |3.4.4| ) but not to the transverse multipole 



-Eo+ which is entirely sensitive to the axial form factor (radius). 

3.4- QCD analysis of charged pion electroproduction 

3.4-1 ■ Effective field theory of the Standard Model In this section, we briefly discuss 
the effective field theory of the Standard Model (SM), chiral perturbation theory, 
without giving any technical details. For that, we refer the reader to a number 
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of reviews 173, |79|, BO, l81|, p2[. At energies and momenta much below the scale 



set by the intermediate vector boson masses, the weak interactions are frozen and 
the SU(3)cxSU(2)LxU(l)y Standard Model simplifies to SU(3)cxU(l), i.e QCD in 
the presence of QED. These two theories behave very differently at low energies. 
QED is characterized by a small coupling constant, oem = 1/137, that means the 
photons and leptons almost decouple and very precise perturbative calculations can 
be performed. On the other hand, the QCD coupling constant becomes large at 
low energies, as = 5's/4vr ~ 1 and the underlying fields (quarks and gluons) are 
confined within hadrons. Consequently, to study the strong interactions, one essentially 
has to consider QCD in the presence of external sources embodying the electroweak 
interactions. Apart from lattice gauge theory attempts, which are plagued by their own 
conceptual problems, no first principle calculations for this theory are possible. However, 
one can make use of the symmetries of QCD and their realizations to explore in a model- 
independent way the chiral dynamics of QCD. This is based on the observation that the 
six quark fiavors fall into two distinct classes, the light (u,d,s) and the heavy (c,b,t) ones 
(light/heavy compared to the typical hadronic scale of 1 GeV). Therefore, to a good 
first approximation one can consider an ideal world with 

T^u = T^d = ms = , rric = nif, = rrit = CO , (37) 

i.e. only the light fiavors are active and the heavy ones decouple. This is called the 
chiral limit of QCD. The corresponding Lagrangian, >Cqq£,, contains only one parameter, 
gs (or, by dimensional transmutation, Aqcd)- -^qcd is highly symmetric because the 
gluon interactions with the quark triplet q = {u, d, s) are fiavor-blind. In particular, 
one can decompose the quark fields in left- and right-handed components. These can 
be independently transformed under SU(3)j^ and SU(3)/j leaving Cqqy) invariant. By 
Noether's theorem, this leads to eight conserved vector, V^ = g7^(A"/2)g, and eight 
conserved axial-vector, A"^ = g7^75(A"/2)(3', currents. From these, one can construct 
conserved charges, 

[Ql.KcB] = [Qt,KcB] (a=l,...,8). (38) 

This is the so-called chiral symmetry of the strong interactions, which has indeed 
been observed long before the advent of QCD and led to current algebra and soft 
pion techniques (as used in the previous subsection). However, this symmetry is not 
necessarily shared by the ground state or the particle spectrum (hidden or spontaneously 
broken symmetry). So what is the fate of the chiral symmetry in QCD? For any vector- 
like gauge theory in the absence of ^-terms, one can show that the vacuum is invariant 



under vector transformations |^3[. Concerning the axial transformations, nature selects 
the Nambu-Goldstone alternative, i.e the ground state |0) is not symmetric and the 
left- and right-handed worlds communicate, as indicated by the non-vanishing quark 
condensate. Because the axial charges do not annihilate the vacuum, the spectrum must 
contain eight pseudoscalar Goldstone bosons, one for each broken generator. These 
particles have zero energy and three-momentum. As the energy /momentum transfer 
decreases, the interaction of the Goldstone bosons with themselves and matter becomes 
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weak, paving the way for a perturbative treatment in energies/momenta or a derivative 
expansion, despite the fact that the strong couphng constant is large. In the real world, 
the 8 lightest hadrons are indeed pseudoscalars (the three pions, four kaons and the eta). 
These particles are light but not exactly massless simply because the quark masses are 
not exactly zero but small. The QCD quark mass term therefore leads to explicit chiral 
symmetry breaking, which can be treated perturbatively. The important observation 
now is that the consequences of the broken chiral symmetry can be explored by means 
of an effective field theory, called chiral perturbation theory (CHPT). In essence, one 
maps the QCD Lagrangian on an effective Lagrangian formulated in terms of the 
asymptotically observed fields, the Goldstone bosons and matter fields (like nucleons), 

>CQCD[g, g, g] -^ C,^[U,d^U,...,M.---,N] , (39) 

where U collects the Goldstone bosons, the quark mass matrix M. parameterizes the 
explicit chiral symmetry breaking and matter fields N can also be included. For a 
deeper discussion of the foundations of CHPT, we refer to the seminal papers by 
Weinberg |^, Gasser and Leutwyler |8^ and Leutwyler [^]. Observables are now 



expanded in powers of momenta (energies), generically called g, and quark masses, at a 
fixed ratio M./q^. Stated differently, one has an underlying power counting |p^ which 



allows one to organize any matrix-element M in terms of tree and loop diagrams for a 
given chiral dimension, symbolically 

M = q''f(^^^,g^ , (40) 

where g is a small momentum or meson mass, / a function of order one, /i a regularization 
scale and gi a collection of coupling constants. Because of chiral symmetry, the counting 
index (chiral dimension) v is bounded from below. One can show that n-loop graphs 
are suppressed by powers of g^*^, thus the leading contribution to any given process 
stems from tree graphs with the lowest order insertions (which is equivalent to current 
algebra). To a given order, the corresponding effective Lagrangian must contain all terms 
compatible with chiral, gauge and discrete symmetries. Beyond leading order, such local 
interactions are accompanied by unknown coupling constants (the gi in equation (0)), 
also called low-energy constants (LECs). While these can in principle be calculated from 
QCD, gi = (yfj(AQCD,'"^c,^6,^t), in practice they must be determined by a fit to some 
data or estimated using some model. Let us illustrate this for the pion-nucleon system 
chirally coupled to external sources like e.g. electroweak gauge bosons. A complete 
one-loop calculation is based on the effective Lagrangian 

Ces = C^^ + C^^ + C^^ + C^^ , (41) 

from which one calculates tree graphs with insertions from all terms >C^^' ' and loop 
diagrams with at most one insertion from the dimension two operators collected in C\j^- 
The lowest order effective Lagrangian contains gA and m (note that the nucleon mass 
might be transformed into higher order 1/m corrections ||88|). At order two, three and 
four one has 7, 23 and 118 LECs, respectively. While these numbers appear large. 
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it should be stressed that many of the dimension four operators contribute only to 
very exotic processes, like three or four pion production induced by photons or pions. 
Moreover, for a given process, it often happens that some of these operators appear in 
certain linear combinations and others simply amount to a quark mass renormalization 
of the corresponding dimension two operator. For example, in the case of elastic pion- 
nucleon scattering, one has 4, 5, and 4 independent LECs (or combinations) thereof at 
second, third and fourth order, respectively, which is a much smaller number than the 
total number of independent terms. Consequently, the often cited folklore that CHPT 
becomes useless beyond a certain order because the number of LECs increases drastically 
is not really correct. It is also important to stress that the values of the dimension 
two (and of some dimension three) LECs can be understood in terms of (low-lying) t- 
channel meson resonance and s-channel baryon resonance excitations |^9| (the so-called 



resonance saturation). For further details, we refer the reader to references [p3, 87 . 



3.4-2. Low-energy theorems for pion electroproduction We now wish to apply the 
machinery of baryon CHPT to the case of pion electroproduction and derive low-energy 
theorems (LETs) . Before doing that, a few clarifying remarks about the meaning of 
such LETs are in order. Following the seminal work of Nambu et al PB| , P^ , pion 
electroproduction LETs have been discussed in the sixties, often using by now outdated 
methods, see [^ |^, ^ |^. As already stressed, only using the methods of CHPT one 
can formulate model-independent statements beyond leading order, provided one is able 
to fix all LECs from data. The modern meaning of LETs is discussed in broad detail in 
reference [0. To third order in the chiral expansion, the various isospin components 
of the transverse and the longitudinal S-wave multipoles can be written as functions of 
the two small parameters /x and v, cf. equation (|3T|) , as (for details, see ||74|| ) 

^i:^(/^, ^) = ^ (-2/^ + /^'(3 + «:.) - Kl + '^.) + rSsi(-'^/^"')| + 0(?') . (42) 

LW(;.,z.) = EW(;.,z.) + ^(/.^-z.){K. + ^S,(-z.;.-^)} + 0(g^), (43) 

4%^ ^) = Ir^ {-2/^ + A^ + f^s) - Kl + '^s)} + Oiq') , (44) 

Jzvrm 

4°^(/i, u) = <)(/., u) + ^ {^? -u)n, + 0{q') , (45) 

iZTXin 



8nm I V 4 8 6 / Sn^r^ 
+ O(g^) , (46) 

with Ky = Kp — Kn = 3.71 the isovector and k^ = Kp + /?„ = —0.12 isoscalar anomalous 
magnetic moment of the nucleon, (/jrAf the strong pion-nucleon coupling constant (which 
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can be used instead of qa by virtue of the Goldberger-Treiman relation) and F^{k'^) is 
vector (charge) form factor of the pion. The functions Sj(— z//yU^), (j = 1,2,3,4) are 
given by 

-.^ P , (2 + P)' 

H2(P) 
S3(P) 



+ 



1 + p ■ 2(l + p)3/2 

2-p (2 + p)2-2 



arccos ■ 
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1 + 
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X 
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These functions are shown in figure ^ for < p < 10. They exhibit a very smooth 
behavior. The influence of isospin breaking, i.e. the charged to neutral pion mass 
difference in the loops, has been discussed in reference WM, here it suffices to say that 
the p-dependence is only very little affected by such effects. Note that the LETs given 
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Figure 6. The functions ^i{p) (i — 1,2,3,4) as defined in equations (p 



above do not contain the full electromagnetic form factors of the nucleon. This is due 
to the power counting of CHPT, from which one concludes that at third order one is 
only sensitive to the normalization of the magnetic form factors, i.e. the respective 
anomalous magnetic moments. At fourth order, one expects sensitivity to the electric 
charge radii of the proton and the neutron. Of particular interest are the LETs for the 
(— ) amplitudes because of the connection to the nucleon axial form factor, as will be 
discussed in the next paragraph. Here, we just mention that the constant C appearing 
in the expression for -£^0+ ^^^ t)e expressed in terms of parameters that can be fixed in 
pion photoproduction, see ||7^. Its explicit form is, however, not needed in what follows. 
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3.4-3. Nucleon axial radius From the LETs just given, we can now derive the 
relation between the electroproduction amphtude and the nucleon axial radius, following 
reference [0. To separate the axial radius, one should consider the slope of the 
transverse multipole. For doing that, one has to expand the function S3(— z//i~^) in 
powers of k'^ = z/m^ and pick up all terms proportional to z/. However, to make the 
argument even clearer, we give one intermediate step. In fact, the direct diagrammatic 
evaluation of the electric dipole amplitude to third order in the chiral expansion includes 
the one-loop diagrams shown in figure |^. These lead to an additional contribution to 
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Figure 7. One-loop diagrams that lead to the axial radius correction. Crossed 
partners are not shown. The solid, dashed and wiggly lines denote nucleons, pions 
and photons, in order. 
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+ 0(g3) . (49) 



Note that in a third order calculation, the electric dipole amplitude is only given 
to second order in small momenta. The constant C in ( PPD subsumes various fc^- 
independent loop corrections together with contact terms and kinematical corrections 
of order 1/m?. In the chiral limit yU = 0, one recovers of course the LET of Nambu et 
al [^ ^ , because the last term of equation (|49|) vanishes identically at zero pion mass 
in the physical region k"^ < 0, compare also equation ( ^61) and the discussion below. 
Matters are, however, different for non-vanishing pion mass, /x 7^ 0. Expanding the 
In-term in (|49D in powers of fc^, the coefficient of the fc^-term in equation (|36|) receives 
an extra contribution, 



i + ir<'-'^> + l!;(''" + ^) + 



e 



-S) 



128F,n^ ttV ^^°^ 

The last term in this equation constitutes a model-independent contribution at order 
k"^ not modified by higher loop or contact terms, which are suppressed by powers of 
/i. Formally, it appears because one cannot interchange the limits M^^ -^ (chiral 
limit) and k"^ -^ (photon point). Therefore, what was believed to be the axial radius 
determined in pion electroproduction was nothing but the modified radius 



(rl) = (rl) + 



64F2 



1 






(51) 
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The numerical value of this additional term is 

A(r^) = (f2)- (r^) = -0.0456 W , (52) 

which translates into a shift of the axial mass of 

AMa = 0.055 GeV , (53) 

which is in agreement with the empirical value given in equation ([TTl) , AM a = 
(0.043 ± 0.026) GeV. This agreement can be considered as one of the premier successes 
of baryon chiral perturbation theory. 

Since the lowest (third) order result for the axial radius discrepancy is quite small, one 
should investigate higher order corrections in /i to find out how reliable the leading order 



result is. This has been reported in |9o] and we essentially follow the arguments given 
there. To order g^, which gives the pertinent next-to-next-to-leading order corrections, 
one has to consider one loop graphs with exactly one insertion from the dimension two 
ttN Lagrangian and counterterms. The LECs related to these will be estimated from 
resonance exchange contributions, here p and A exchanges. Since we are only interested 
in the transverse part of the transition matrix element for 7*p -^ ir^n at threshold, we 
introduce an auxiliary quantity E{k'^), that is proportional to the transverse threshold 
S-wave multipole, 

T e = ^ a . e E{e) = 47rz(l + p) a ■ e Ef^^l^, . (54) 

In the chiral limit, the corresponding current algebra result becomes exact and gives, 
when expanded in k"^ 

E{k') = 1 + y (rl) - ^ (^„ +\^ + O(fc^) , (55) 

with 

(rl) = {r\) + 0{p') , L= n^ - 4^ + 0(/^') ' (56) 

the nucleon axial radius and the neutron magnetic moment in the chiral limit. Here the 

o 

superscript '°' denotes quantities in the chiral limit, Q =Q [1 + 0(^M,d)]. Note that the 
axial radius (like the isocalar electromagnetic radius) is analytic in the quark masses. 
We now want to calculate all tree and loop graphs up to order q^ which contribute to 
the slope E'{Q) = dE{0)/dk'^ and are proportional to /i° and/or /i. The quantity -^'(O) 
is the sum of the desired squared axial radius and a host of other terms. Among these 
are contributions from the Born graphs including electromagnetic form factors. Since 
these are already contained in the standard analysis of the pion electroproduction data, 
such terms need not be considered further. The discrepancy subsumes all loop and 
counterterm effects which go beyond the form factors. After some lengthy calculation, 
one arrives at 
,~2\ /^2\ 3 f^ 12^^ , 3M^ , 3c+(7r-4)M^ 



(^a) (rA) g4^2 (1 ^2 ) + 64^ p2 + 327rF2 



-i^('"^-^-S4)-->'--^'°'' (-' 
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The first term in ( |57| ) is the leading order result already given before. The combination 
of LECs, c+ = —8ci+4c2 + 4c3 — (7^/2777,, can be related to the isospin-even ttN scattering 
length aj^_ via c^ = 87rF^a^_^/M^ (which is correct to the order we are working here). 
The last two terms in (|57| ) represent the counterterm contributions at order g^, which 
have been identified with p(770) and A(1232) exchange contributions, 

3(1 + Kp)M^ 



6e;(o) : 



IGTT'^gAfnF^ 



\/2vn?'m\ 



rriA — m 
-2m{Y + Z + 2YZ)-2mAiY + Z + AYZ) , (5^ 



with Kp ~ 6 the tensor-to-vector ratio of the pNN couplings and k* = (71 ~ 5 is 
the ^AN coupling constant. The so-called off-shell parameters Y and Z (which are 
nothing but higher order LECs in an effective field theory including the delta) have 



been constrained in reference [^ from the delta contribution to the proton's magnetic 
polarizability, Y = —0.12, and from the 033 ttN scattering volume, Z = —0.3. Putting 
all pieces together, one finds for the axial radius discrepancy, 

{r\) - {r\) = (-4.6 + 0.5) ■ 10"^ fm^ . (59) 

Note that there are large individual terms in the loop and the counterterm contributions 
which cancel each other. Also, this correction is fairly sensitive to the badly known 
isoscalar ttN scattering length, so it is fair to assign an uncertainty to this correction 
which is as large as it central value given in equation (|59|) . However, the only meaningful 
quantity is the total sum of all terms of order q'^, which can also be considered verified 
by the MAMI experiment |^^. Taken that experiment face value, one could turn the 
argument around and give a bound on aj^_. We refrain from doing that here but 
stress again the intricate interplay of seemingly unrelated processes, linked by the chiral 
structure of QCD. It should also be noted that the momentum transfers in the MAMI 
experiment were too high to allow for a direct and safe extraction of i?o+ 5 therefore 
the data were analyzed in the framework of an effective Lagrangian model with the 
electromagnetic nucleon form factors, the pion vector form factor and the axial nucleon 



form factor at the appropriate vertices [^, ^ . A further measurement at lower energy 



and lower momentum transfer should therefore be studied for feasibility. 

3.4-4- Pion charge radius The charge (vector) radius of the pion is a fundamental 
quantity in hadron physics. It can essentially be determined in two ways. One method 
is pion scattering off electrons (or electron-positron annihilation into pion pairs), this 
leads a pion root-mean-square (rms) radius of [^ 

{rDlf^ = (0.663 ± 0.006) fm , (60) 

if one insists on the correct normalization of the pion charge (vector) form factor, 
F^{0) = 1 (in units of the elementary charge e). A more recent determination of the 
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pion vector radius from low-momentum space- and time-like form factor data based on 



a very precise two-loop chiral perturbation theory representation ||100|| gives (two-loop 
contributions to the pion form factors were first considered in reference [|101|| ) 

(r2)j/2 = (0.661 ± 0.012) fm , (61) 

consistent with the value given above. Electromagnetic corrections considered in 
reference [102| reduce this value by about one percent. This number can be understood 



semi-quant it at ively in a naive vector meson dominance picture, (r^)y = ^/6/Mp ~ 
0.63 fm. The second method is based on charged pion electroproduction, 7*p — > n^n. 
The unpolarized cross section in parallel kinematics decomposes into a transversal and 
a longitudinal piece, as detailed in section |3.4.2| . While the former is sensitive to the the 



nucleon axial radius, the latter is quite sensitive to the pion form factor, i.e. to the pion 
radius for small momentum transfer. A recent measurement at the Mainz Microtron 
MAMI-II led to a pion radius of [||] 

(r2);/2 = (0.74 ± 0.03) fm, (62) 

which is a sizeably larger value than the one obtained from vre scattering. It was 
hinted in reference p^ that their larger value for the pion radius might be due to 



the inevitable model-dependence based on the Born term approach to extract the pion 
radius. It was also stated that there might be an additional correction obtainable form 
chiral perturbation theory as it is the case for the nucleon axial radius discussed before. 
Indeed, there is such a similar kind of correction for the pion radius [ |105|| . This new term 
modifies the momentum dependence of the longitudinal S-wave amplitude Lq, and leads 



one to expect an even larger pion charge radius than the one given in equation 
It is conceivable that higher order corrections yet to be calculated or contributions 
from higher multipoles will completely resolve the discrepancy between the pion radius 
determined from vre scattering on one side and from charged pion electroproduction on 
the other. Let us be more specific now. As shown in equation (0), the pion vector form 
factor Fj^(fc^) appears in the expression for the longitudinal multipole. The pion form 
factor has the following low-energy expansion, 

F^{e) = l+^-{rl)ye + 0{k'). (63) 

Consequently, to separate the pion radius, one should consider the slope of the 
longitudinal multipole. For doing that, one has to expand the function S4(— z/yU~^) 
in powers of k"^ = vrn^ and pick up all terms proportional to v. This gives: 



dL[-^ 



(9fc2 



^9.N / 1 l_ , i3_ 1 2\ ^ I 1 /^16 



fc2=o 327rmlM2 mM^ Sm^ ?,^ ^' 32F^\7i^ 



(-^)^ + ^^b-l)+^(^^)}- (64) 



The first four terms are standard | |106| ], they comprise the conventional dependence 



on the pion vector radius, recoil effects and the dominant chiral limit behavior of the 
slope of the longitudinal multipole. The strong 1/M^ chiral singularity stems from 
the fc^-derivative of the pion-pole term which appears already at leading order in the 
chiral expansion (this is the leading contribution from the induced pseudoscalar form 
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factor). The last term in equation (0) originates from the so-called triangle and 
tadpole (with three pions coupling to the nucleon at one point) diagrams which are 
known to play a prominent role in pion photo- and electroproduction (see figure |^. 
The formal reason for the appearance of this new, model-independent contribution at 
order k"^ is that one cannot interchange the order of taking the derivative at fc^ = 
and the chiral limit M^^ —>■ 0. Consequently, all determinations of the pion radius 
from electroproduction (based on tree-level amplitudes including nucleon and pion form 
factors) have "measured" the modified radius, 

(av = (av + ^(i|-l). (65) 

The novel term on the right-hand-side of equation (p5| ) amounts to 0.266 fm^, a bit 
more than half of the squared pion rms radius, {r'^)v — 0.44 fm^. Therefore, from the 
longitudinal multipole alone, one expects to find a larger pion radius if one analyses 
pion electroproduction based on Born terms, 

{7l)y = (0.44 + 0.26) fm^ = (0.83 fm)^ , (66) 

which is even larger than the result of the Mainz analysis, cf. equation (jB^). We 



point out, however, that the contribution of the pion radius to the derivative of the 
longitudinal multipole is a factor of ten smaller than the one from the first three terms 



in the curly brackets in equation (|6^) . Therefore, a fourth order analysis is certainly 
needed to further quantify the "pion radius discrepancy" . Furthermore, the pion form 
factor contribution to the longitudinal cross section is also present in higher multipoles. 
In fact, it is known that the convergence of the multipole series for the pion pole term 
is slow. One should therefore also investigate such effects for these higher multipoles 
or directly compare the predictions of complete one-loop calculation with the data of 
the longitudinal electroproduction cross section. For the purpose of demonstrating the 
significance of chiral loop effects the /c^-slope is, however, best suited. What has been 
shown here is that as in the case of the nucleon axial mean square radius, the pion loops, 
which are a unique consequence of the chiral symmetry of QCD, modify the naive Born 
term analysis and should be taken into account. 

4. Induced pseudoscalar form factor and coupling constant 

4-1- QCD analysis of the induced pseudoscalar form factor 

Here, we show that one can give an accurate prediction for the induced pseudoscalar 
form factor Gp{t) and coupling constant gp in terms of well-known physical parameters, 
following reference ||107|| (for a similar analysis, see [|108|| ). For doing that, one exploits 



the chiral Ward identity of two-fiavor QCD, 



d^ (^T^Tsyg) =rhqi-i-,T''q (67) 
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with m the average hght quark mass (note that isospin-breaking effects can be safely 
neglected here). Sandwiching equation (p7| ) between nucleon states, one obtains ||109| 



mGAit) + -^Gp{t) = 2mBmh ^j^^ (68) 

Am M^ — t 

The pion pole in equation (|68D originates from the direct coupling of the pseudoscalar 



density to the pion, (0|gi75r"g|7r ) = 5 G^ |85|. The residue at the pion pole t = M^ is 



|T09| 



rh G^ g^N = 9-kN F^ M^ (69) 

with Qttn the strong pion-nucleon coupling constant. The chiral expansion of the axial 
form factor and of the function h{t) to order q^ takes the form 

GA{t)=gA(l+^-^t] , h{t) = const - ^t (70) 

V 6 y gA 

with dis a low-energy constant (we use the notation of reference ||110|| ) and one does 
not need to specify the constant since it is not needed explicitly in the following. The 
reason for the linear dependence in equation (^) is the following. The corresponding 
form factors Ga(^) and h{t) have a cut starting at t = (SM^r)^ which in the chiral 
expansion first shows up at two-loop order 0{q^). Therefore, the contribution to order 
g^ must be polynomial in t. Furthermore, from chiral counting it follows that the t^- 
terms are related to order q^ of the full matrix-elements. Putting pieces together, we 
arrive at 

mgA + mgA^^t+-—Gp{t) = -— — -t + g^NF^ + (71) 

6 4m M^ — t gA 

using 2mB gAm= M^{g^NF^ + 0{M^)). At t = 0, equation (|7l]) reduces to the 
Goldberger-Treiman discrepancy (GTD) |p,09| , |111 | 



9Am = g^^ fJi + '^ M^) , (72) 

which also clarifies the meaning of the low-energy constant rfig. We will return to a 
more detailed discussion of the GTD in section ^. Finally, Gp{t) can be isolated from 
equation ([7T|), 

Gp{t) = ^^^ -\gAm- {r\) + 0(t, M^J . (73) 

A few remarks are in order. First, notice that only physical and well-determined 
parameters enter in this formula. Second, while the first term on the right-hand-side 
is of order g~^, the second one is 0{q^) and the corrections not calculated are of order 
g^. Third, the first term is the celebrated pion pole term, which stems from the direct 
coupling of the pion field to the pseudoscalar source. It clearly dominates this form 
factor for momentum transfers of the order of a few pion masses. Also, because of the 
pion pole, the structure of this form factor is not of the common form, i.e. its low- 
energy expansion is not a polynom characterized by some normalization and a radius. 
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Because of the Ward identities, it is linked to the axial form factor as witnessed by the 
appearance of the axial radius in equation (^). For the pseudoscalar coupling gp, this 
leads to 

Qp = — :r-^ ?r Qa m,, m ir a) . 74 

^ Ml + 0.88M2 3 y^ /^ \ A/ y ! 



Indeed, this relation has been derived long time ago by Adler and Dothan ||112|| and 



by Wolfenstein ||113[|. Wolfenstein used a once-subtracted dispersion relation for the 



right-hand-side of equation (p5D (weak PCAC). It is gratifying that the ADW result 
can be firmly based on the systematic chiral expansion of low energy QCD Green 
functions. In chiral perturbation theory, one could in principle calculate the corrections 
to equation ([7^) by performing a two-loop calculation while in Wolfenstein's method 
these could only be estimated. To stress it again, the main ingredient to arrive at 
equation ([7^) in CHPT is the linear t-dependence in equation ([701). Since we are 
interested here in a very small momentum transfer, like e.g. for the case of ordinary muon 
capture, t = — 0.88M^ ~ — 0.5M^, curvature terms of order t^ have to be negligible. 
If one uses for example the dipole parameterization for the axial form factor, the t^- 
term amounts to a 1.3% correction to the one linear in t. The masses m, M^ and 
M^ = Mn-+ are accurately known and so are F^ and qa 0- The situation concerning 
the strong pion-nucleon coupling constant is less favorable. The methodologically best 



determination based on dispersion theory gave g'^j^/Air = 14.28 ±0.36 [p.l4]| , more recent 



determinations seem to favor smaller values ||115|| . We use here (/^at = 13. 10 ±0.35. The 



value for the axial radius has been discussed in section |2.1| , we use here the number 
obtained from neutrino scattering. Putting pieces together, we arrive at 

gp = (8.74 ± 0.23) - (0.48 ± 0.02) = 8.26 ± 0.16 . (75) 

The uncertainties stem from the range of g.,rN and from the one for {r\) for the first and 
second term, in order. For the final result on gp, we have added these uncertainties in 
quadrature. A measurement with a 2% accuracy of gp, as intended by the upcoming 
experiment ^^, could therefore cleanly separate between the pion pole contribution and 
the improved CHPT result. This would mean a significant progress in our understanding 
of this fundamental low-energy parameter since the presently available determinations 



have too large error bars to disentangle these values, see section p.2| . In fact, one 
might turn the argument around and eventually use a precise determination of gp to 
get an additional determination of the strong pion-nucleon coupling constant which 
has been at the center of much controversy over the last years. To summarize, we have 
shown that the chiral Ward identities allow to predict the induced pseudoscalar coupling 
constant entirely in terms of well-determined physical parameters within a few percent 
accuracy. As already noted by Wolfenstein |p,13|| , an accurate empirical determination 
of this quantity therefore poses a stringent test on our understanding of the underlying 
dynamics which is believed to be realized in the effective low-energy field theory of 
QCD. 
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4-2. Chiral expansion of ordinary muon capture 



As detailed in section |2.2| , measurements of OMC on liquid hydrogen allow one to extract 
a value for the induced pseudoscalar coupling constant that agrees with theoretical 
expectations. The theoretical description of the hadronic physics involved is fairly 
transparent, we follow reference ||116|| and refer to that paper for many details. However, 



as stressed in particular by Ando, Myhrer and Kubodera ||117|| , there are some difficulties 
with the accepted picture of translating the atomic rates into the one for liquid hydrogen. 
This topic will be taken up at the end of this section. First, we consider the OMC 
(atomic) decay rates calculated to third order in the chiral expansion and also to third 
order in the so-called "small scale expansion", which is an effective field theory with 



explicit A(1232) degrees of freedom (for a detailed introduction and review, see ||118|| ) 



In the Fermi approximation of a static W~ field, i.e. the gauge boson propagator is 
reduced to a point interaction (since the typical momenta involved are much smaller 
than the W mass), 

M,-,^.^n = -MO^c = {u^\w;\fx)tl^ [(n|Klp) - {n\A-\p)] , (76) 

where the vector current matrix element {n\V~\p) (the so-called vector correlator) is 
parameterized in terms of the isovector vector Dirac and Pauli form factors F^2j which 
are known empirically (see e.g. ||119| , |120|| ) and have also been studied in baryon 



CHPT ||111| , |121| , |122|| . The axial correlator contains the desired dependence on the 
induced pseudoscalar form factor. Introducing the Fermi constant Gp, the square of the 
spin- averaged invariant matrix element is defined via 

1 ^|_^OMC|2^1^^^^2^^(a)^.._ (77) 



4 



crcr'SS' 



with Vud = 0.974 the pertinent CKM matrix element and the explicit forms of the 
symmetric leptonic L^,^, see also equation (|1^), and the hadronic tensor H'^'^ are given 



in |116| . Approximating now the muonic atom wave-function dependence by its value 
at the origin (thus neglecting the small third order effects due to the form factors), the 
Is Bohr wave function of the muonic atom reads 

HO)u = -^ , (78) 



with the reduced mass fi = {mM^-)/{m + M^~) and a the fine-structure constant. We 
can therefore calculate the spin-averaged rate of ordinary muon capture via 



OMC 



with J,y, Jn appropriate normalization factors. Evaluating this to third order, O(g^) 
one obtains 



OMC = ( ^VAJ -bi.tD -a.» + 0(l/m^) j X s"^ = 181.7 x s" 




(80) 
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which shows a very fast convergence. Retaining the delta only reduces the small third 
order correction by 5%. The reason for the nice stability of perturbative calculations 
for OMC in the physical world of small finite quark masses is of course the fact that 
contributions of order n are suppressed by (Mj/A^)""" , with i = 7r,/i and A^ ~ m ~ 
1 GeV. This spin-averaged OMC scenario is mostly of theoretical interest. In nature the 
weak interactions show a very strong spin-dependence, which leads to quite different 
decay rates depending on whether the captured Is muon forms a singlet or a triplet spin- 
state with the proton [|123|] , where the singlet is the usual state (1/a/2)(| Tj i) ^ I i; T)) 
in terms of the muon and proton spins and the triplet accordingly, see also figure |^ and 
the discussion in section |2]^. Although the hyperfine splitting between the two levels 
"only" amounts to 0.04 eV, the occupation numbers of the levels due to thermal and 
collision induced processes tend to be far from statistical equilibrium. In order to make 
any contact with experiment, the singlet /triplet rates need to be calculated separately. 
For the total capture rates of singlet and the triplet states in the muonic atom, one find 
the following decomposition into leading, next-to-leading and next-to-next-to-leading 
order pieces 

T'^^^ = (957 - 245 GeV/m + (30.4 GeV/m^ - 43.17) + 0{l/nV^)) x s"^ 

= 687.4 X s"^ , 
r*(5JP(. = (10.3 + 4.72 GeV/m - (1.22 GeV/m^ + 1.00) + 0{l/m^)) x s"^ 

= 12.9 X s~^ , (81) 

displaying the dramatic spin-dependence due to the V-A structure of the weak 
interactions in the Standard Model. These numbers correspond to a value of q-^nn = 
13.05 and include the pion pole corrections, cf. equation (^). Since Gp contributes 
negatively to the singlet rate a larger value of gnNN leads to a smaller value for the rate: 
-'^OMC ~ 681.9 X s~^ for Qt^nn = 13.4. Similarly, neglecting the pion pole corrections 
leads to Fq^^ = 676.1 x s^^. In Equation (^) we have split the third order term (third 
and fourth terms in parenthesis) into the contribution from the 1/m^ corrections and 
the terms stemming from the various radii which lead to the Q^-dependence of the form 
factors. It is very interesting to note that these two contributions more or less cancels 
themselves in the case of the singlet term. It is thus extremely important to perform 
a consistent chiral expansion. The existing theoretical predictions for these rates are 
collected in table S. For comparison, in a relativistic Born model one obtains the simple 



Table 2. Calculated atomic singlet and triplet rates in s from BHM |11£] 



AMK [ |117| , Primakoff [|2J and Opat |123(| . Note that the latter two have used smaller 



values for qa and that AMK used Qt^n — 13.4. 

BHM BHM AMK AMK Primakoff Opat 
NLO NNLO NLO NNLO 

r«ing 711 687.4 695 722 664 ±20 634 

r*"P 14.0 12.9 12.2 11.9 11.9 ±0.7 13.3 
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nice formula for T^^q ||125|| : 



"'^s -■ {Q.2ZQF^{ql) + 0.5513 ^2^(50") + lQAAGA{ql) ' 0.2834 Gp(go'))' 



OMC 

-1 



~ 683 X S-' (82) 

where q^ = — 0.88M^ is the momentum transfer corresponding to muon capture at rest 
(again, a rescahng of this formula to the present day value of qa has been performed). 
This formula though very appealing should not be used, being in contradiction with 
the modern viewpoint of power counting. The good agreement with the CHPT result 
is purely accidental. The problem arises from the fact that Tq^q is a rather sensitive 



quantity as can be seen for example in equation (82). Indeed the terms proportional to 
-^2'('?o) ^^^ Gp{lo) ^^6 of the same order of magnitude but have different signs so they 
have a tendency to cancel each other rendering the values of r'^"^^ rather sensitive to the 
exact values of these two quantities. So far we have only considered OMC for the case of 
muonic atoms. For the case of a liquid hydrogen target one also has to take into account 
the possibility of muon capture in a muonic molecule p/xp, which can be formed via the 
reaction p/i + pep -^ pfip + e + 124 eV. In such a molecule the muon can be found in the 
ortho (O) (spin of the protons parallel) or para (P) (spin of the protons antiparallel) 
spin state relative to its two accompanying protons. The decay rates of these molecular 
states can be calculated from the singlet /triplet rates of the muonic atom via 

Tp = 27p-(3rtrip + Tsing) , (83) 

To = Pl/2 ri/2 + P3/2 r3/2 , (84) 

with 7p (70) denoting the ratio of the probability of finding the negative muon at 
the point occupied by a proton in the para-muonic (ortho-muonic) molecule and the 
probability of finding the negative muon at the origin in the muonic atom. The 



wavefunction corrections are taken to be 70 = 0.500, 7p = 0.5733 ||126|| . We note 
that the para molecular state is often referred to as the statistical mixture, as it 
corresponds to the naively expected occupation numbers of the muonic atom. For a 
precise calculation of the ortho molecular state on the other hand one first has to know 
the exact probabilities pi/2,P3/2 for the muonic molecule being in a total spin S=(l- 
l/2)=l/2 or a total spin S=(l+l/2)=3/2 state, with pi/2 > 0.5 |^, see also figure |. 
The corresponding decay rates are given by |124, TB 



ri/2 — 270 ( TFsing + Trtrip 

13/2 = 270 Ftrip (85) 

Theoretical calculations of the spin-effects in the muonic molecule [^ suggest pi/2 ~ 
1, P3/2 ~ which leads to the values: 



o^c = 208 X s-i 

OMC _ Ann _ _ _ cm X. .-l 



r^^^^ = 493---519xs-V (86) 

where the range given in r^^'^' corresponds to 0.95 < pi/2 < 1. We have allowed here 
for a 5% uncertainty in the occupation numbers to show the sensitivity of our results on 
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this quantity. This point was first made in reference ||117|| . Since Fsing ^ Ftrip for OMC, 
-pOMC lyfj-Yia, out to be roughly proportional to pi/2- We note that our number for capture 
from the molecular ortho state agrees very well with the most recent measurement 
V^Q^ = (531 ± 33) X s~^ |5^. We stress again that a direct use of the liquid transition 
formula a la reference [Q is only meaningful if one accounts for the time structure of 
the beam in the experiment, see also 



4-3. Chiral expansion of radiative muon capture 

The pioneering TRIUMF RMC result, cf. equation ([151) , spurred a lot of theoretical 
activity. While radiative muon capture had already been calculated in phenomenological 
tree level models a long time ago, see e.g. [Q, |123| , |127| , |128| , |129|| , heavy baryon chiral 
perturbation theory was also used at tree level including dimension two operators ||130 



and to one loop order ||131|| . The resulting photon spectra are not very different from 
the ones obtained in the phenomenological models, the most striking feature being 
the smallness of the chiral loops and polynomial third order corrections ||131|| , hinting 
towards a good convergence of the chiral expansion. At present, the puzzling result 
from the TRIUMF experiment has not fully been explained, but it appears now that 
the discrepancy does not come only from the strong interactions but rather is also 
related to the distribution of the various spin states of the muonic atoms, as detailed 
below [|116| , |117|] . It is therefore mandatory to sharpen the theoretical predictions for the 
strong as well as the non-strong physics entering the experimental analysis. Therefore, 
RMC was also analyzed in the framework of the small scale expansion [|118|| , which allows 
to systematically include the A resonance into the effective field theory. Although Ando 



and Min ||131|| have already shown that the RMC process possesses a well behaved chiral 
expansion up to N^LO, it has been noted quite early ||132|| that one should reanalyze 
RMC in a chiral effective field theory with explicit A degrees of freedom. This is 
due to the fact that the A-resonance lies quite close to the nucleon and therefore, in a 
delta-free theory like baryon CHPT, could lead to unnaturally large higher order contact 
interactions which would spoil the seemingly good chiral convergence. Stating this in the 
language of (naive) dimensional analysis, it suggests the possibility of corrections of the 
order of 30% due to the small nucleon-delta mass splitting, M^/{m^ — m) r^ 3M^/m. 
Such a study including explicit delta degrees of freedom was performed in ||116|| . It 
should also be mentioned that phenomenological models have claimed for a long time 
that the A contribution does not exceed 8% in the photon spectrum for photon energies 
above 60 MeV ||129| . We will now proceed to discuss general features of the chiral and 
the small scale expansion for RMC and then critically reexamine the way the TRIUMF 
experiment was analyzed. 



4-3.1. General results In the static approximation for the W-boson, the pertinent 
matrix element for RMC decomposes into two terms. 



M 



M P^i^iin^ 



{u,\W^\fi)tjj^ [{n\TV ■ e*V-\p) - {n\TV ■ e*A-\p)] 



M'w 
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87) 



where the first one contains the vector-vector (VV) and vector-axial (VA) correlator. 
The hadronic part of the second term is, of course, identical to the OMC matrix- 
element. Following ||116|| , we give these to second order in the chiral and the small scale 
expansion. The third order terms have been worked out by Ando and Min ||131|| and 
found to be small. Working in the Coulomb gauge for the photon and making use of 
the transversality condition e* ■ A; = 0, one finds (the pertinent momenta were already 
given in equation (0)) 



{n\TV-e*V-\p) 



(2) 
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with UJ = V ■ k, S^ the covariant nucleon spin-vector, v^ the nucleons four-velocity 
(see e.g. reference ||111|| for a more detailed discussion of the underlying heavy 
baryon formalism) and i? = 1 in QCD. Furthermore, Kg is the nucleon isoscalar 
anomalous magnetic moment, A = m^ — rn the delta-nucleon mass splitting, g-j^NA 
and hi are the leading vrA^A and 7A^A coupling constants. Numerically, Qt^n/^ x bi = 
1.05 X 12 = 12.6 ||118|| . The SU(2) gauge coupling g2 is related to the Fermi constant 
via Gf = 5'2 V^/(8M^). Note that the vector-vector correlator is free of delta effects at 
next-to-leading order, i.e. the leading A(1232) effect only appears in the vector-axial 
correlator. The latter terms constitute the difference between the chiral and the small 
scale expansions for RMC evaluated to NLO. The factor R multiplying the Born term 
contributions proportional to the induced pseudoscalar form factor has been introduced 
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for the later discussion. Note again that these expressions do not contain the full form 
factor Gp{t) but only its leading part (the CHPT correction given in equation (^) does 
not appear at this order). One can easily check from the continuity equations satisfied 
by the correlators that gauge invariance is satisfied in the above equations. We are now 
in the position to calculate the decay rates. Consider first the spin-averaged case. The 
square of the matrix-element equation (|57| ) can be written as a sum of four terms, with 
both photons coming either from the hadronic or the leptonic side and two mixed terms, 
i.e. 

e'GlVl 



,(90) 



aa'ss'W L \ AA' 

with A, A' the photon helicities. Explicit expressions for the various tensors are not given 
here because they are lengthy and not illuminating. The total decay rate is given by: 

[^(o)!,^ r . r— / /m^-^(i-cos0) 



167r4 Jo Jo ° V V m 

aa'ss'W 

with uj = kf) the photon energy (note that in this and the following subsection we 
exceptionally use the symbol u for the photon energy to facilitate the comparison with 
the literature). The direction of the photon defines the z-direction and 6 in equation (0) 
is the polar angle of the outgoing lepton with respect to this direction. The maximal 
photon energy is given by 

...^M„(:.^)(:.^)-'. (.2) 

Furthermore, the energy of the outgoing lepton follows from energy conservation, 

m m 

First we discuss the (academic) spin-averaged RMC scenario, which allows for a 
comparison with previous calculations: 

^fZ^^.. = (66.0 + 18.7 + C(l/m^)) x 10"^ s"^ = 84.7 x 10"^ s"^ (CHPT) 

TsTnav. = (66.0 + 20.4 + C(l/m2)) x lO^^^ s"^ = 86.4 x lO^^^ s"^ (SSE). (94) 

Both the CHPT and the SSE results suggest a good convergence for the chiral expansion, 
as expected from dimensional analysis. Note that the leading order capture rates in both 
calculations are identical, as A(1232) related effects only start at sub-leading order. Our 



leading order result also agrees with the calculation of reference ||130|| . For the case of 
muonic atoms we obtain the following decay rates in the singlet /triplet channel (note 
the reversal of the relative size of the singlet to triplet contribution as compared to the 
OMC case), 

Tf^^ = (12.7 - 18.7 GeV/m + 0{l/m^)) x lO^^ s"^ = 3.10 x lO^^ s"^ (CHPT) 
^fmg^ = (12.7- 18.3 GeV/m + C(l/m2)) x lO^^^s"^ = 2.90 x lO^^^s"^ (SSE) 
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rf^ljc ^ (119- 3.86 GeV/m + 0(l/m2)) x 10"^ s"^ = 112 x IQ-^s"^ (CHPT) 

Tf^^ = (119 - 1.80 GeV/m + 0{l/m'^)) x 10"^ s^^ = 114 x 10"^ s"^ (SSE) (95) 

It should be stressed that for the total numbers given the kinematical factors were 
not expanded in powers of 1/m since in case of the small singlet, the contribution 
from the terms starting at order 1/m? can not be neglected. Next, we address the 
complications for RMC due to the presence of muonic molecules in the liquid hydrogen 
target. According to equation (p3D, we can easily determine the capture rate from the 
molecular para state 

r^^^c _ gg 2 [86.4] X 10-^s-i (CHPT [SSE]) . (96) 

Consider now the molecular ortho state, which turns out to dominate in the recent RMC 
experiment from TRIUMF p2| , |53[ . One obtains for pi/2 = 1: 

r^MC _ 3Q_4 pg g] ^ 10-3 s-1 (CHPT [SSE]) , (97) 

i.e in both cases the delta effects are fairly small. Due to the triplet dominance in 
RMC (as opposed to the singlet dominance in OMC) Pq^*-' is now roughly proportional 
to (1 — 3/4^1/2) which leads to a big sensitivity of the RMC capture rate to the exact 
occupation numbers of the relative molecular sub-states. For example, a 5% uncertainty 
in the occupation numbers pi/2 = 0.95, P3/2 = 0.05 would lead to a 13% change in the 
ortho capture rate P^^*^ ~ 35 x 10~^s~^. We will discuss the implications of this 
uncertainty when we compare our results with the measured photon spectrum from 
TRIUMF in the next subsection. For the total capture rate in the TRIUMF experiment 

J- RMC - /5 J- sing + Jo J- O + JP ^ P \d^) 

with fs = 0.061, fo = 0.854, fp = 0.085 [§g, §3] one would obtain P^^^^^*^^ = 
(34.3 [34.8] + 0(l/m2)) x lO'^s'^ in HBCHPT [SSE]. This leads to a relative branching 
ratio Q^ = PRMc/roMC 

Q^ = 6.8 [6.9] X 10"^ + 0(l/m2) (CHPT [SSE] ) . (99) 

Unfortunately the full relative branching ratio is not accessible in experiment, as one 
has to use a severe cut on the photon energies due to strong backgrounds. In the 
TRIUMF experiment only photons with an energy u; > 60 MeV were detected. We 
therefore now move on to a discussion of the photon spectrum dT/du. It can be obtained 
straightforwardly from equation (0). The calculated photon spectra to second order in 
CHPT and SSE come out very similar. In figure |^ we show the second order SSE results 
for the singlet, triplet, para and ortho states. The relative difference between the second 
order CHPT and SSE calculation for all states are of the order of a few percent, showing 
explicitly the small role of A (1232) in RMC. These results are very similar to the ones 
found by Beder and Fearing [|6^, |129|| for the spectra and the relative contribution from 
the spin-3/2 resonance, although their calculation is based on a very different approach. 
Even the result for the small singlet is comparable though not identical to the one of 
Beder and Fearing. We conclude that the strong interaction theory for RMC is under 
control. 
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Figure 8. Photon spectra for RMC for the singlet, triplet, para (statistical) and ortho 
states of the /i — p system calculated to second order in the small scale expansion. 



4.3.2. Discussion of the TRIUMF result The photon spectra just discussed allow in 
principle to determine the induced pseudoscalar form factor. The TRIUMF result for 
gp is obtained by multiplying the terms proportional to the pseudoscalar form factor 
with a constant denoted R (the momentum dependence is assumed to be entirely given 
by the pion pole). The value of R is then extracted using the model of Fearing et al 
to match the partial rate for photon energies larger than 60 MeV. If we perform such a 
procedure, we get a similar shift in the partial photon spectra. It is, however, obvious 
from the analysis presented above that such a procedure is not legitimate. By artificially 
enhancing the contribution ~ gp (to simulate this procedure, we have introduced the 
factor R in equation (|89D and similarly one must multiply {n\A~\p) by i?), one mocks 
up a whole class of new contact and other terms not present in the Born term model. To 
demonstrate these points in a more quantitative fashion, we show in figure ^ the partial 
branching fraction in comparison to the one with gp enhanced by a factor 1.5 and a 
third curve, which is obtained by increasing gp only by 15% and slightly modifying 
some parameters related to the delta contribution. This is shown by the dashed line 
and it shows that such a combination of small effects can explain most (but not all) 
of the shift in the spectrum. This is further sharpened by using now the neutral pion 
mass of 134.97 MeV instead of the charged pion mass, leading to the dotted curve in 
figure 1^. Since the pion mass difference is almost entirely of electromagnetic origin, one 
might speculate that isospin-breaking effects should not be neglected. Furthermore, as 
discussed above a slight change in the occupation numbers pi/2 and P3/2 would also lead 
to an increase in the ortho capture rate which could close the gap between the empirical 
and theoretical results. For example the dashed curve in figure ^ would be moved from 
0.42 to 0.48 while the dotted one would go from 0.50 to 0.54 with pi/2 = 0.95 and 
P3/2 = 0.05 at tu = 60 MeV. The situation is reminiscent of the sigma term analysis, 
where many small effects combine to give the sizeable difference between the sigma term 
at zero momentum transfer and at the Cheng-Dashen point. Another point against the 
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Figure 9. Photon spectra for RMC for the branching ratios of the singlet, ortho and 
para states as used in the TRIUMF analysis. Solid line: Prediction of the small scale 
expansion to order e^. Dash-dotted line: Same as the solid line but with gp scaled by 
a factor 1.5. Dashed line: Various small modifications as explained in the text. Dotted 
line: Same as the dashed line but using the neutral instead of the charged pion mass. 



simple rescaling comes from OMC Indeed if this rescaling holds for RMC it should also 
hold for OMC. We thus have performed a similar calculation in OMC. Taking the same 
value for R, one would obtain 



'- OMC 



172.8 X s-\ 

OMC, R=1.5 



-\sing, R=1.5 
OMC 



634.6 X s 



-1 



itrip, R=1.5 
OMC 



18.9 X S 



-1 



(100) 



leading to Fq"^""' "^~^-^ = 477 x s~^, which is lower than the error bars of the experimental 
result from reference [^. As expected, the singlet and triplet capture rates are much 
more sensitive to the details of the interaction than the total rate. As a conclusion 
we note that the effect of enhancing the capture rates in RMC via setting R = 1.5 
leads to a strong reduction of the corresponding OMC rates leading to conflicts with 
the experimentally determined ortho capture rate. To summarize this discussion, we 
have pointed out that two effects in particular have to be investigated in more detail: 
(1) the occupation numbers of the atomic structure in muonic atoms/molecules need to 
be carefully re-examined, and (2) the N^LO calculation should be redone including all 
isospin breaking effects because of the sensitivity to the exact pion mass in the pion-pole 
contributions, for example. The sum of these small effects should explain the observed 
photon spectrum, as we believe that the proper hadronic/weak physics part is well 
under control by now. A simple rescaling of the pseudoscalar coupling constant should 
no longer be considered. 



5. Octet Goldberger-Treiman discrepancies 

In section ^ we already encountered the Goldberger-Treiman discrepancy (GTD), i.e. 
the deviation from the Goldberger-Treiman relation (GTR), which is exact in the chiral 
limit. This deviation is parameterized in terms of the quantity A^^ given by 



A. 



mgA 



F^g^NNiM^] 



g^NN^M^ 



(101) 
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using the GTR mgAiS^) = FT^gnNNiO). Inserting the PDG values for qa, -PV and 
m = [iTLp + m„)/2, one obtains 

A^ = 0.014 (0.040) for q^nn = 13.05 (13.40) . (102) 

From naive dimensional analysis one can can estimate A^^ ~ (M^r/A^)^ = (Mt^/AttFj^)^ ~ 
0.015, so this favors the smaller pion-nucleon coupling constant. Note, however, that if 
one chooses the p-meson mass as the relevant hadronic scale, A^ = Mp, this estimate 
increases to {M^^/MpY — 0.033. One can further sharpen this argument by extending 
the analysis to flavor SU(3) and also investigate consistency of the strange to light 
quark mass ratio as obtained from the Goldstone boson spectrum. This is based on the 
observation that while in SU(2) the GTD is given in terms of one LEG, in chiral SU(3) 
this quantity is given in terms of a sum of quark masses times an SU(3) octet operator. 
This leads to a relation between various GTDs in the octet, the so-called Dashen- 
Weinstein relation ||133|| . This can be explored more systematically in the framework of 
CHPT. We follow here essentially reference ||134]| and the recent update given in ||135 



(see also ||136|| ). In the baryon octet, one calculates matrix-elements of the axial-currents 
u-^tTl^d and ifj^^s at zero momentum transfer, 

{B'{p)\A^\B{p)) = UB' [(7f'^7M75 + --]uB, (103) 

and the pertinent axial coupling constants g^ ^ can be measured in semi-leptonic 
hyperon decays (/5-decays), B -^ B'iu^. Extending the analysis of section^ to the SU(3) 
case, one can derive a variety of GTR's and the chiral corrections to these. We focus 
here on the three cases for which sufficient empirical information for the corresponding 
strong meson-baryon coupling constants is available, i.e. np, Ap and E~n transitions. 
To linear order in the quark masses, the chiral Ward identity relating the divergence of 
the axial current to the pseudoscalar density sandwiched between baryon states leads 
to 

(m„ + mp)gj = 2F^g^NN + (m„ + md)iJ"P(0) + 0(mj) , 

(mA + mp)gf = -v^F^(7a + (m„ + m,)if^P(0) + 0(m;) , 

(ms- + mO^F" = 2FKgj, + (m„ + m,)iJ^""(0) + 0(mj) , (104) 

with gA = —gApK- and gs = g-E-nK- /V^. Furthermore, the functions H^^' (t) are SU(3) 
generalizations of h(t) in equation (pBD, given by 

{p\ui-f5d\n) = [H''^P(t) + 7r - pole] Upij^Un , 

(p|m75s|A) = [H^'\t) + K- pole] VTs^a , 

{n\ui'j5s\T.~) = H^'''{t) + K -pole Unh^uj:- . (105) 

Note that we are using the convention for the axial charges so that g^ is positive. The 
H functions admit a chiral expansion, 

H^P{0) = F + D + O{m,) , 

/7Ap(0) = -yi(F+|Z})+O(m,), 

H^~''{0) = ~F + D + O{mg) . (106) 
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Here, F = 0.46 and D = 0.80 are SU(3) axial-vector couplings. Therefore we have 
three relations with two unknowns, from that we can deduce a sum rule in terms of the 
GTD's (A^, A^, A|-) and the quark masses, sometimes called the Dashen-Weinstein 
relation (DWR), 

A. = ^ -^^ (Vs^-^^Ai + ^^^A| ] , (107) 

h^ m + rus V 9-kNN Q-kMN J 

in terms of the kaon discrepancies Aj^' , 



^~' 2FKgKKN{Ml^ 



Ai = 1 - 






A^ = l- ^ "^^ "^ '\Zn,' " ■ (108) 



(m „ + ms-)(-(7r" (0)) 
2FKgKMMl) 

and rh = (m.^ + mii)/2 the average light quark mass. There are two ways to 
explore the DWR. First, consider the quark mass ratio appearing in equation ( |107| ) 
known, 2m/{m + m^) = M'^/Mj^ ~ 0.08. Note that higher order strong as well as 
electromagnetic corrections do not change this leading order estimate appreciably, see 
e.g. ||137||. Then one can ask the question whether three GTDs are mutually consistent. 



that is whether the DWR is fulfilled. Another way of exploring this relation is to 
determine the three discrepancies from data and then use the DWR to deduce the quark 
mass ratio in /{in + nis). As input, we use the PDG values for the masses and axial 
coupling constants together with gxANiM].) = 13.7 ±0.4 and gKJ^NiM]^) = 3.9 ±0.4 as 



obtained from the analysis of hyperon-antihyperon production at LEAR ||138|| . These 
numbers are close to the expectations from flavor SU(3), gxAN = 13.5 and gxi^N = 4.3, 
respectively, and in rough agreement with the analysis of ^N scattering data given in 
reference ||139|| . This leads to the following values of the kaon discrepancies 

A^ = 0.16 ± 0.06 , A| = 0.18 ±0.11 . (109) 

These numbers conform to the expectations from SU(3) symmetry breaking, Ak — 
(M^/A^)^ = 0.18, but the theoretical uncertainties might be larger as the ones given, 
see the discussion at the end of this section. Using now the quark mass ratio as deduced 
from the Goldstone boson masses and the central values of the GTDs, one obtains 
A^ = 0.016, which favors the low value of the pion-nucleon coupling constant. If, on 
the other hand, one uses the empirical values of the GTDs, see equations (|102yi09| ), one 
obtains for the quark mass ratio 

m^^i 22.4 ± 6.8 for g^,,^ = 13.05 

m I 7.4 ± 3.3 for ^^^^v = 13.40 ' ^ ^ 

i.e. the smaller value of g^NN is consistent with the standard scenario of chiral symmetry 
breaking, that is a large value of the quark condensate, B = (0|gg|0)/F^ ~ 1.5 GeV. 
One could be content with this consistency since the so determined quark mass ratio 
is in good agreement with the one obtained from analyzing the Goldstone boson mass 
spectrum based on the large condensate scenario. However, as already stressed in ||134||, a 
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better and more reliable determination of the strong coupling constants Qkan and gKY:N 
is called for. In this context it is worth to stress that extraction of these couplings from 
kaon electroproduction data usually leads to very different values, see e.g. [|140|| and the 



discussion in reference ||141|| . Note, however, that most of these kaon electroproduction 



analyses are based on simple Born terms (including resonance excitations and form 
factors) and are therefore plagued with large theoretical uncertainties. In the light of 
the new and coming data on kaon electropoduction off protons and deuterium from 
ELSA and CEBAF, it would be very valuable to apply the theoretical framework of 
references [|142|, |143|] , which combines chiral perturbation theory, unitarity and coupled 



channel dynamics (see also [|144| , |145|] ). The presently existing uncertainties for these 



coupling constants inhibit a precise test of the chiral symmetry breaking pattern of 
QCD. 

6. Summary and outlook 

The nucleon as probed with the weak axial current can be parameterized in terms of two 
form factors, the axial, GA{t), and the induced pseudoscalar, Gp{t), one. In this short 
review, we have shown that precise theoretical methods based on the symmetries of QCD 
exist for extracting these fundamental observables from experiment. The axial form 
factor can be well described by a dipole, Gyi(t) = (1 — t/M^)~^. The dipole mass Ma 
can be translated into an axial root-mean-square radius of (?"i)^^^ = 0.67±0.01 fm. This 
value is consistently obtained from (anti)neutrino scattering off protons (or light nuclei) 
and charged pion electroproduction off protons. Clearly, more precise electroproduction 
data in the threshold region would be welcome to further pin down this quantity. The 
induced pseudoscalar form factor is dominated by the pion pole, but the small corrections 
to this leading order result have been calculated. Existing data from ordinary muon 
capture are consistent with these theoretical expectations but have too large error bars 
to cleanly test the chiral dynamics of QCD. We have argued that the result of the 
pioneering TRIUMF radiative muon capture experiment should be taken cum grano salis 
due to some assumptions in the analysis that are inconsistent with the power counting 
underlying the effective field theory of the Standard Model. However, it is fair to say that 
more theoretical as well as experimental effort is needed for drawing a final conclusion. 
The momentum-dependence of the induced pseudoscalar form factor is dominated by 
the pion pole which has only been tested in one electroproduction experiment so far. 
Also, a better determination of strong kaon-hyperon-nucleon coupling constants Qkan 
and Qk'sn would allow for an additional stringent bound on the light to strange quark 
mass ratio, based on the derivations from the octet Goldberger-Treiman relations. All 
this shows that precision experiments in hadron and nuclear physics indeed help to 
unravel the mysteries of QCD at energies where one really has to deal with strong 
interactions. Therefore, more pion and kaon electroproduction as well as muon capture 
experiments are called for. 



Axial structure of the nucleon 37 

Acknowledgments 

We are grateful to Norbert Kaiser and Thomas Hemmert for fruitful collaborations, 
Fred Myhrer for sharing his insight into the theory of muon capture and Simon Sirca for 
detailed explanation of the MAMI data and analysis. Useful communications from Barry 
Holstein, Jacques Martino and Claude Petitjean are acknowledged. We thank Josef 
Speth for initiating this project. We also thank the Deutsche Forschungsgemeinschaft 
and the ELSA/MAMI Schwerpunktprogramm for funding the workshop on "Hadron 
Form Factors" in Bad Honnef in April 2001, which triggered this write-up. This work 
was supported in part by the NSF, grant 34201. 

References 



[1 

[2: 

[3 

[4; 

[6 
[T. 

[8 
[9 

[lo; 

[11 

[12 
[13 
[14 
[15 
[16 

[17; 
[is: 

[19 



[20 

[21 
[22 
[23 
[24 
[25 
[26; 
[27; 
[28 
[29 
[30 
[31 
[32 
[33; 



Ulf-G. Meifiner, Phys. Rcpt. 161 (1988) 213. 

S. Weinberg, Phys. Rev. 112 (1958) 1375. 

M. Morita, R. Morita and K. Koshigiri, Nucl. Phys. A577 (1994) 387c. 

D.E. Groom et al (Particle Data Group), Eur. Phys. J. C5 (2000) 1. 

J. Reich et al , Nucl. Inst. Meth. A440 (2000) 535. 

H. Abele, Nucl. Inst. Meth. A440 (2000) 499. 

I.S. Towner and J.C. Hardy, in Symmetries and Fundam,ental Interactions in Nuclei^ E.M. Henley 

and W.C. Haxton (eds.) (World Scientific, Singapore, 1995), pp. 183-249. 
G. Fanourakis et al , Phys. Rev. D21 (1980) 562. 
L. A. Ahrens et al , Phys. Rev. D35 (1987) 785. 
L. A. Ahrens et al , Phys. Lett. B202 (1988) 284. 
S. J. Barish et al , Phys. Rev. D16 (1977) 3103. 
K. L. Miller et al , Phys. Rev. D26 (1982) 537. 
W. A. Mann et al , Phys. Rev. Lett. 31 (1973) 844. 
N. J. Baker et al , Phys. Rev. D23 (1981) 2499. 
T. Kitagaki et al , Phys. Rev. D28 (1983) 436. 
T. Kitagaki et al , Phys. Rev. D42 (1990) 1331. 
M. Holder et al , Nuovo Cim. ALVII (1968) 338. 
R. L. Kustom et al , Phys. Rev. Lett. 22 (1969) 1014. 

D. Perkins, in: Proceedings of the 16"^ International Conference on High Energy Physics, 
J. D. Jackson, A. Roberts (eds.). National Accelerator Laboratory, Batavia, Illinois, 1973, 
Vol. IV, 189. 

A. Orkin-Lecourtois and C. A. Piketty, Nuovo Cim. AL (1967) 927. 
S. Bonetti et al , Nuovo Cim. A38 (1977) 260. 
N. Armenise et al , Nucl. Phys. B152 (1979) 365. 
I. Budagov et al , Lett. Nuovo Cim. II (1969) 689. 
A. Liesenfcld et al , Phys. Lett. B468 (1999) 20. 

E. Amaldi et al , Nuovo Cimento A65 (1970) 377. 
E. Amaldi et al , Phys. Lett. B41 (1972) 216. 
E.D. Bloom et al , Phys. Rev. Lett. 30 (1973) 1186. 
P. Brauel et al , Phys. Lett. B45 (1973) 389. 

A. del Guerra et al , Nucl. Phys. B99 (1975) 253. 

A. del Guerra et al , Nucl. Phys. B107 (1976) 65. 

P. Joos et al , Phys. Lett. B62 (1976) 230. 

A. S. Esaulov, A. M. Pilipenko, Yu. I. Titov, Nucl. Phys. B136 (1978) 511. 

M. G. Olsson, E. T. Osypowski and E. H. Monsay, Phys. Rev. D17 (1978) 2938. 



Axial structure of the nucleon 38 



[34] S. Choi, Etude experimentale et theorique de Velectroproduction de pions pres du seuil et mesure 
des facteurs de forme axial et pseudo-scalaire, Saclay report DAPNIA/SPhN 93 18, 1993. 

D. Baumann, U. Miiller et al (Al Collaboration), forthcoming. 
Y. Nambu and D. Lurie, Phys. Rev. 125 (1962) 1429. 
Y. Nambu and E. Shrauner, Phys. Rev. 128 (1962) 862. 
Y. Nambu and M. Yoshimura, Phys. Rev. Lett. 24 (1970) 25. 
G. Furlan, N. Paver and C. Verzegnassi, Nuovo Cim. A LXX (1970) 247. 

C. Verzegnassi, Springer Tracts in Modern Physics 59 (1971) 154. 

G. Furlan, N. Paver and C. Verzegnassi, Springer Tracts in Modern Physics 62 (1972) 118. 
N. Dombey and B. J. Read, Nucl. Phys. B60 (1973) 65. 

B. J. Read, Nucl. Phys. B74 (1974) 482. 

G. Benfatto, F. Nicolo and G. C. Rossi, Nucl. Phys. B50 (1972) 205. 
G. Benfatto, F. Nicolo and G. C. Rossi, Nuovo Cim. A14 (1973) 425. 
S. Weinberg, Phys. Rev. Lett. 4 (1960) 575. 
E.J. Bleser et al , Phys. Rev. Lett. 8 (1962) 288. 
J.E. Rothberg et al , Phys. Rev. 132 (1963) 2664. 

E. Bertolini et al , Proc. Intern. Conf. on High-energy physics (CERN, Geneva 1962) p. 421; 
S. Focardi et al , Proc. Intern. Conf. on Fundamental aspects of weak interactions (BNL, Upton, 
1964), p. 278. 

V.M. Bystritskii et al , Sov. Phys. JETP 39 (1974) 19. 

G. Bardin et al , Nucl. Phys. A352 (1981) 365. 

G. Bardin et al , Phys. Lett. 104B (1981) 320. 

D.D. Bakalov, M.P. Faifman, L.I. Ponomarev and S.I. Vinitsky, Nucl. Phys. A384 (1982) 302. 

S.-I. Ando, K. Kubodera and F. Myhrer, in preparation. 

D.V. Bahn et al , PSLproposal R-97-05. 

D.F. Measday, The Nuclear Physics of Muon Capture^ Phys. Rept., in print. 

Ch. Briancon et al , Nucl. Phys. A671 (2000) 647. 

P. Ackerbauer, Phys. Lett. B417 (1998) 224. 

G.H. Miller, Phys. Lett. 41B (1972) 50. 

F.R. Kane, Phys. Lett. 45B (1973) 292. 

P. Guichon, Phys. Rev. C19 (1979) 987. 

G. Jonkmans et al , Phys. Rev. Lett. 77 (1996) 4512. 

D.H. Wright et al , Phys. Rev. C57 (1998) 373. 

D. Beder and H.W. Fearing, Phys. Rev. D39 (1989) 3493. 
S. Choi et al , Phys. Rev. Lett. 71 (1993) 3927. 
P. Bartsch et al (Al Collaboration), MAMLproposal Al/1-98. 

C. Llewellyn Smith, Phys. Rept. 3C (1972) 261. 

E. Fischbach et al , Phys. Rev. D15 (1977) 97. 
J. Kim et al , Rev. Mod. Phys. 53 (1981) 211. 
J. Collins, F. Wilczck and A. Zee, Phys. Rev. D18 (1978) 242. 
D.B. Kaplan and A.V. Manohar, Nucl. Phys. B310 (1988) 527. 



[35 

[36; 

[37 
[38 
[39; 
[40' 

[41 
[42 
[43; 
[44; 
[45 

[46; 

[47 
[48; 
[49; 



[50 
[51 
[52 
[53 
[54 
[55 
[56 
[57; 
[58 
[59 
[6O; 
[61 
[62 
[63 
[64 
[65 
[66 
[67 



[69; 

[70 
[71 
[72 
[73 
[74 
[75 
[76 
[77 
[78 
[79 
[80 
[81 



W.M. Alberico, S.M. Bilenky and C. Maieron, subm. to Phys. Rept., |iep-ph/0102269 

E. Amaldi, S. Fubini and G. Furlan, Springer Tracts in Modern Physics 83 (1979) 6. 

V. Bernard, N. Kaiser, T.-S.H. Lee and Ulf-G. Meil3ner, Phys. Rept. 246 (1994) 315. 

F.A. Berends, A. Donnachie and D.L. Weaver, Nucl. Phys. B4 (1967) 1. 

J. Bernstein, S. Fubini, M. Cell-Mann and W. Thirring, Nuovo Cimento 17 (1960) 757 

S. Coleman, A.spects of Symmetry (Cambridge University Press, Cambridge, 1985). 

H. Georgi, Ann. Rev. Nucl. Part. Sci. 43 (1993) 209. 

Ulf-G. Meifiner, Rep. Prog. Phys. 56 (1993) 903. 

A. Pich, Rep. Prog. Phys. 58 (1995) 563. 

G. Ecker, Prog. Part. Nucl. Phys. 35 (1995) 1. 



Axial structure of the nucleon 39 



[82 
[83; 
[84 

[85 
[86 

[87; 



[89; 

[90' 
[91 

[92 
[93; 

[94; 
[95 

[96; 

[97; 



100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 

114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 



V. Bernard, N. Kaiser and Ulf-G. MciBncr, Int. J. Mod. Phys. E4 (1995) 193. 

C. Vafa and E. Witten, Nucl. Phys. B234 (1984) 173. 
S. Weinberg, Physica 96A (1979) 327. 

J. Gasser and H. Leutwyler, Ann. Phys. 158 (1984) 124. 

H. Leutwyler, Ann. Phys. 235 (1994) 165. 

Ulf-G. Meifiner, in: Boris loffe Festschrift At the Frontier of Particle Physics - Handbook of QCD, 

Vol. 1, pp. 417-501, Ed. M. Shifman (World Scientific, Singapore, 2001). 
E. Jenkins and A.V. Manohar, Phys. Lett. B255 (1991) 558. 
V. Bernard, N. Kaiser and Ulf-G. MeiBner, Nucl. Phys. A615 (1997) 483. 
S. Fubini, G. Furlan and C. Rosetti, Nuovo Cim. 40 (1965) 1171. 
Riazuddin and B.W. Lee, Phys. Rev. 146 (1966) B1202. 
S.L. Adler and F.J. Gilman, Phys. Rev. 152 (1966) B1460. 
S.L. Adler, Ann. Phys. (NY) 50 (1968) 189. 

G. Ecker and Ulf-G. Meifiner, Comm. Nucl. Part. Phys. 21 (1995) 347. 
V. Bernard, N. Kaiser and Ulf-G. Meifiner, Phys. Rev. Lett. 69 (1992) 1877. 

A. Schmidt, Pion-Elektroproduktion nahe der Schwelle in chiraler Storungstheorie, PhD thesis, 
TU Miinchen (1995). 

K. Blomquist et al (Al Collaboration), Z. Phys. A353 (1996) 415. 

D. Drechsel, L. Tiator, J. Phys. G: Nucl. Part. Phys. 18 (1992) 449. 

S. Amendolia et al (NA7 Collaboration), Nucl. Phys. B277, 168 (1986). 
J. Bijnens, G. Colangelo and P. Talavera, JHEP 9805 (1998) 014. 
J. Gasser and Ulf-G. Meifiner, Nucl. Phys. B357 (1991) 90. 

B. Kubis and Ulf-G. Meifiner, Nucl. Phys. A671 (2000) 332. 

V. Bernard, N. Kaiser, J. Kambor and Ulf-G. Meifiner, Nucl. Phys. B388 (1992) 315. 

T. Becher and H. Leutwyler, Eur. J. Phys. C9 (1999) 643. 

V. Bernard, N. Kaiser and Ulf-G. Meifiner, Phys. Rev. C62 (2000) 028201. 

S. Scherer and J.H. Koch, Nucl. Phys. A534 (1991) 461 and references therein. 

V. Bernard, N. Kaiser and Ulf-G. Meifiner, Phys. Rev. D50 (1994) 6899. 

H.W. Fearing, R. Lewis, N. Mobed and S. Scherer, Phys. Rev. D56 (1997) 1783. 

J. Gasser, M.E. Sainio and A. Svarc, Nucl. Phys. B307 (1988) 779. 

N. Fettes, Ulf-G. Meifiner and S. Steininger, Nucl. Phys. A640 (1998) 199. 

V.Bernard, N. Kaiser, J. Kambor and Ulf-G. Meifiner, Nucl. Phys. B388 (1992) 315. 

S.L. Adler and Y. Dothan, Phys. Rev. 151 (166) 1267. 

L. Wolfenstein, in: High-Energy Physics and Nuclear Structure, ed. S. Devons (Plenum, New 

York, 1970) p.661. 
G. Hohler, in Landolt- Bornstein, vol. 9b2, ed. H. Schopper (Springer, Berlin, 1983). 
R.A. Arndt, LI. Strakovsky and R.L. Workman, Phys. Rev. C52 (1995) 2246. 
V. Bernard, T.R. Hemmert and Ulf-G. Meifiner, Nucl. Phys. A686 (2001) 290. 
S.-I. Ando, F. Myhrer and K. Kubodera, Phys. Rev. C63 (2000) 015203. 
T.R. Hemmert, B.R. Holstein and J. Kambor, J. Phys. G24 (1998) 1831. 
P. MergeU, Ulf-G. Meifiner and D. Drechsel, Nucl. Phys. A597 (1995) 367. 
H.-W. Hammer, Ulf-G. Meifiner and D. Drechsel, Phys. Lett. B385 (1996) 343. 
V. Bernard, H.W. Fearing, T.R. Hemmert and Ulf-G. Meifiner, Nucl. Phys. A635 (1998) 121. 
B. Kubis and Ulf-G. Meifiner, Nucl. Phys. A679 (2001) 698. 
G.I. Opat, Phys. Rev. 134 (1964) B428. 
H. Primakoff, Rev. Mod. Phys. 31 (1959) 802. 
A. Santisteban and R. Pascual, Nucl. Phys. A260 (1960) 392. 
W.R. Wessel and P. Phillipson, Phys. Rev. Lett. 13 (1964) 23. 
D. Beder, Nucl. Phys. A258 (1976) 447. 
H.W. Fearing, Phys. Rev. C21 (1980) 1951. 
D. Beder and H.W. Fearing, Phys. Rev. D35 (1987) 2130. 



Axial structure of the nucleon 40 



[130; 
[131 
[132 
[133 
[134 
[135 
[136 
[137; 
[138 

[139 
[140; 
[141 
[142 
[143 
[144 
[145 



T. Meissner, F. Myhrer and K. Kubodera, Phys. Lett. B416 (1998) 36. 

S.-I. Ando and D.-P. Min, Phys. Lett. B417 (1998) 177. 

N.C. Mukhopadhyay, in Proceedings of WEIN-98 (Santa Fe, NM, 1998), ^iucl-th/9810039 

R. Dashen and M. Weinstein, Phys. Rev. 188 (1969) 2330. 

N.H. Fuchs, H. Sazdjian and J. Stern, Phys. Lett. B238 (1990) 380. 

J.L. Goity, R. Lewis, M. Schvellinger and L.-Z. Zhang, Phys. Lett. B454 (1999) 115. 

B.R. Holstein, Few Body Syst. Suppl. 11 (1999) 116. 

H. Leutwyler, Phys. Lett. B378 (1996) 313. 

R.G.E. Timmermans, T.A. Rijken and J.J. de Swart, Phys. Lett. B257 (1991) 227; Nucl. Phys. 

A585 (1995) 143c. 
A.D. Martin, Nucl. Phys. B179 (1981) 33. 
R.A. Adelseck and B. Saghai, Phys. Rev. C42 (1990) 108. 

H. Haberzettl, C. Bennhold, T. Mart and T. Feuster, Phys. Rev. C58 (1998) 40. 
Ulf-G. Meii3ner and J.A. Oiler, Nucl. Phys. A673 (2000) 311. 
J.A. Oiler and Ulf-G. Meifiner, Phys. Lett. B500 (2001) 263. 
N. Kaiser, T. Waas and W. Weise, Nucl. Phys. A612 (1997) 297. 
E. Oset et al , Prog. Part. Nucl. Phys. 44 (2000) 213. 



